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Preface 

Accurate  predictions  of  the  signal  to  noise  ratios  that  will 
cause  an  Amplitude  Comparison  Monopulse  Radar  to  break  angle  tracking 
lock,  are  at  present,  not  available.  This  thesis  models  the 
Amplitude  Comparison  Monopulse  Radar  in  terms  of  a  tracking  loop 
using  the  Maximum  Likelihood  Criteria.  From  this  model,  probability 
density  functions  are  desired  for  the  estimated  angle  off  boresight. 

A  criteria  for  predicting  break  lock  is  derived  from  the  probability 
density  functions. 
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Abstract 


An  amplitude  comparison  monopulse  radar  is  modeled  using 
additive  channel  and  system  noise  to  the  received  signals.  The 
amplitude  of  the  incoming  signal  and  the  angle  off  boresight  are 
estimated  under  the  Maximum  Likelihood  Criteria.  An  ensemble  of 
estimates  of  the  angle  off  boresight  are  used  to  derive  probability 
density  functions  for  the  estimate  angle  off  boresight.  From  these 
probability  density  functions,  a  criteria  for  predicting  breaklock 
is  derived. 


ix 


I  INTRODUCTION 


Background 

Accurate  predictions  of  the  signal  to  noise  ratios  that  will 
cause  an  Amplitude  Comparison  Monopulse  Radar  to  break  angle  tracking 
lock,  are  at  present  not  available.  To  enable  accurate  predictions  to 
be  made,  probability  density  functions  for  the  estimated  angle  off 
boresight  need  to  be  derived  for  various  signal  to  noise  ratios. 

Techniques  for  analysis  and  prediction  of  break  lock  for  Phase 
Lock  Loops  have  been  developed  and  well  documented  (References  11  and 
12).  Also,  the  analysis  of  a  Phase  Lock  Loop  system  modeled  under 
Maximum  Likelihood  Estimate  Criteria  for  simple  estimation  problems, 
has  been  documented  (Reference  11) . 

The  Amplitude  Comparison  Monopulse  Radar  has  been  modeled  under 
an  approximation  to  the  Maximum  Likelihood  Criteria  by  Bakut  (Refer¬ 
ence  1) .  Bakut  shows  how  all  system  noise  and  channel  noise  can  be 
lumped  to  one  normally  distributed  noise  source  which  is  additive  to 
the  received  signal. 

Cohen  and  Steinmetz  (Reference  3)  have  modeled  the  antenna 
functions  of  the  Amplitude  Comparison  Monopulse  Radar  as  the  sum  of 
sines  and  cosines  but  this  is  an  approximation  of  the  antenna  pattern 
within  the  3db  beam  width. 

Objective 

To  gain  accurate  probability  density  functions  of  the  angle  off 

A 

boresight  estimate  (8),  an  accurate  model  of  the  radar  needs  to  be 
derived.  The  antenna  functions  are  to  be  modeled  without  any 
approximations.  The  radar  will  be  modeled  under  the  Maximum  Likeli¬ 
hood  Estimate  Criteria.  A  discrimination  model  to  convert  error 
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voltage  to  angle  estimate  off  boreslght  Is  to  be  modeled.  The 
resulting  tracking  loop  probability  density  functions  of  the  angle 
estimate  with  their  relevant  mean  and  variances  are  to  be  plotted. 

Assumptions 

The  scope  of  this  thesis  was  restricted  by  the  following 
assumptions  in  relation  to,  the  noise,  the  radar  platform,  the  radar 
model  and  the  target  movement. 

The  noise  is  assumed  to  be  white  gaussian  distributed  noise.  As 
the  noise  is  additive  to  the  sum  and  difference  channels  of  the  radar 
it  is  also  assumed  to  be  independent  from  pulse  to  pulse  and  from 
channel  to  channel. 

The  model  is  restricted  to  one  dimension  only  and  total  de¬ 
coupling  from  a  second  dimension  is  assumed.  In  deriving  the  model, 
noise  is  assumed  to  be  zero  mean  as  well  as  white  independent 
gaussian  distributed.  The  radar  platform  is  assumed  stationary  and 
consequently  eliminates  modulation  of  the  radar  input  due  to  platform 
movement. 

The  target  is  assumed  to  have  two  separate  movements;  stationary 
and  a  constant  velocity.  When  the  target  is  assumed  to  be  stationary, 
then  this  also  includes  the  case  where  the  radar  is  tracking  a 

target  which  is  moving  slowly  in  relation  to  the  total  number  of  j 

pulses  Integrated .  I 


II  MODEL  OF  AMPLITUDE  COMPARISON 
MONOPULSE  RADAR 

General 

.  Through  the  use  of  two  (minimum  number)  offset  antennas. 

Figure  1,  an  Amplitude  Comparison  Monopulse  Radar  provides  range  and 
angle  information.  The  range  information  is  derived  from  the  signal 
return  delay.  Angle  information  is  derived  by  comparing  the  addition 
and  subtraction  of  the  return  signal's  output  from  the  two  offset 
antenna,  and  then  estimating  the  angular  error  of  the  target  to  the 
antenna  system  boresight  axis .  The  general  block  diagram  of  an 
Amplitude  Comparison  Monopulse  Radar  is  shown  in  Figure  2. 


Figure  1.  Offset  Antenna  Configuration 


OJ t ?) 


Figure  2.  General  Block  Diagram  of  Amplitude  Comparison  Monopulse 
Radar 

Noise 

To  account  for  system  noise  and  channel  noise,  noise  is  lumped 
and  modeled  as  independent  gaussian,  zero  mean  and  added  independently 
to  the  sum  channel  and  difference  channel  to  form  the  X  signals.  The 
input  signals  to  the  estimator,  the  XE  (t)  and  XA(t)  signals  contain 
amplitude,  range  and  angle  off  system  boresight  information  and  are 
modeled  as  follows: 

XT (t)  -  AgE (9)S(t+x^)  +  nE (t) 

XA(t)  -  AgA(6)S(t+T2)  +  nA (t) 

Received  Signal  S(t+T) 

As  this  thesis  is  concerned  primarily  with  angle  tracking  the  t 
(range  estimate)  problem  will  be  ignored.  This  consequently  enables 
the  S(t+r)  portion  of  the  X  signals  to  be  evaluated  as  S(t).  There- 
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fore,  the  X  signals  become: 

XE(t)  -  AgE(9)S(t)  +  nE(t) 

XA(t)  -  AgA(9)S(t)  +  nA(t) 

The  A  component  of  the  received  signal  Is  a  random  process  of  un¬ 
known  density.  The  S(t)  signal  will  be  modeled  in  two  methods. 

(1)  As  a  transmitted  pulse  of  unit  energy  with  the 
received  function  constant 

S(t)  =  u(t)  CKtfT 

0  else 

This  represents  a  return  from  a  stationary  target. 

(2)  As  a  transmitted  pulse  of  unit  energy  with  angular 

movement  and,  consequently,  the  received  antenna 
functions  as  a  function  of  0(t) 

S(t)  =  ^  u(t)  0<t<T 

0  else 

0(t)  =  Pt  +  Q 

Where  P  is  the  rate  of  change  and 
Q  is  the  starting  angle. 

In  both  cases  the  amplitude  of  the  incoming  signal  A  is  the 
random  variable. 

The  typical  pulse  duration  of  a  monopulse  tracking  radar  is 
0. 2p8.  To  aid  in  the  computer  implementation  of  the  radar  model, 
the  pulse  has  been  scaled  to  20ms  with  all  other  timing  scaled 


accordingly 


S(t) ;  0  <  t  £  20ms. 

Antenna  functions  r£(6),  gA(9) 

A  pencil  beam  of  7.5°  for  each  antenna  was  assumed.  To  create 
the  pencil  beam  pattern,  a  square  aperture  with  uniform  aperture  field 
distribution  was  also  assumed.  A  representation  of  the  field  is 
Figure  3. 


A 


K - M — *1 


Ea  ■  Eox/x/<a; (y)<b 
0  else 

Figure  3.  Square  aperture  with  uniform  aperture  field  distribution 
The  fourier  transform  (Reference  4)  of  this  aperture  field 
gives  the  function  ft.  as: 

j  m  x  4abEo  sin  (koa  sin  8  cos  4>)  sin  (kob  sin  6  sin  <1>) 
(koa  sin  0  cos  $)  (kob  sin  6  sin<t>) 

In  the  far  field, 

E(r)  *  1  e  ^or  $  Cos  0sin  <{>(-fx)  +  0  fx  cos  <j> 

/it  r 


.  IT  3-TT  _ 

in  $  -  2*  jf  Plane; 


fx 


4abEo  sin  (kob  sin  6) 
kob  sin  6 


and; 


1koe"jkor 

E(r)  -  fx  +  cos  0  4 

where  <f>  ■ 


for  the  high  gain  antenna  case; 

„  - — >  m  jkoe  ^kor  4abEo  sin  (kob  sin  8) 
r  2trr  kob  sin  0 


Normalizing  for  E(o); 

E(r)  m  sin  [kob  sin  (j 
E(o)  kob  sin  0 

a' 

The  normalized  radiation  pattern  for  each  antenna  is 


Figure  4 


kj|=< 


Figure  4.  Normalized  radiation  pattern  for  each  antenna  3db 
beam  width  is  assumed  as  7.5°. 

Therefore, 

sin  (koa  sin  3.75°) 

.707  =  - r-» - .  ,  7Eo — c 

koa  sin  3.75 


.707  =  sine  (koa  .0654) 

This  gives  koa  *  13.  6874 tt 

koa  =  43.0003097  (2) 

The  two  antennas  form  relative  main  radiation  beams  as  in 


Figure  5.  The  two  antenna  boresight  axes  have  a  separation  of  d. 
For  this  analysis  both  beams  will  be  referenced  to  the  antenna 
system  boresight  axis. 
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Figure  5.  Relative  Antenna  A  and  B  main  radiation  beams 
From  Equation  1 


Antenna  A 

e  e 
e  e 

An 

(0) 

sin 

koa  sin  (0  + 

koa  sin  (  6  +  y) 

Antenna  B 

e  e 

E  0 

£1 

(0) 

m  8  in 

jkoa  sin  (e  -  ^)j 

koa  sin  (0  _ 

If  the  separation  between  the  two  lobes  A  and  B  is  small 
compared  to  the  focal  length,  then  squint  angle  -r  can  be  computed 
approximately  from  the  antenna  parameters  (References  3  and  6)  by 
the  following: 

d  =  Kp  tan  ^  j 

.  *  _  -1  .4 

d  *  • D  tan  »  c , £ 

5.546 
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This  produces  a  radiation  pattern  as  Figure  7. 


Figure  7.  Radiation  pattern  for  gA (8 ) 

Estimator  9 

At  this  point  in  the  Monopulse  Radar,  (output  of  the  I.F. 
section,  input  to  the  estimator) ,  the  information  signals  are  the 
XE(t)  and  XA(t)  signals.  From  these  two  signals  6,  (angle  the  target 
is  off  boresight),  must  be  estimated. 

A  way  of  estimating  0  is  to  maximize  the  a  posteriori  density 
of  p0/f(XE,  XA,  a)  (0/f(XE,  XA,  A)),  (References  10  and  11),  as  in 
Figure  8.  If  the  maximum  is  interior  to  the  allowable  range  of  0 
and  In  p0/f(XE,  XA,  a)  (0/f(X£,  XA,  A))  has  a  continuous  first 
derivative  then  a  condition  for  a  maximum  can  be  obtained  by  differen¬ 
tiating  In  p  0/f(XI,  XA,  a)  (0/f(X£,  XA,  A))  with  respect  to  0  and 
setting. the  result  to  zero. 


Figure  8.  An  a  posteriori  density 
Let  f(XI,  XA,  A)  =  Y 
From  Bayes  Rule  (Reference  Van  Trees) 

p9/y  (em .  nojmLjjs) 

Therefore, 

In  p6/y  (0/Y)  ■  In  y/B  (Y/0)  +  In  6(0)  -  lny  (Y) 

As  the  main  point  of  interest  here  is  to  find  the  value  of  0 
where  the  left  hand  side  is  maximum,  the  last„term  on  the  right  hand 
side  is  not  a  function  of  0  and  can  be  ignored. 

Therefore,  to  maximize  the  a  posteriori  density  In  9/y  (0/Y), 
need  only  maximize; 

lnp  y/0  (Y/0)  +  In  0  (9) 


Nothing  is  known  about  the  a  priori  knowledge  of  0(6)  and  therefore 

if  this  is  assumed  to  be  a  uniform  distribution  within  the  range  of  6 

then  the  a  priori  knowledge  is  a  constant  K. 

3  In  K  -  0 


BOt 


30 
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Therefore,  to  maximize  In  0/y  (0/Y)  is  equivalent  to  maximizing 
In  y/9  (Y/0  )  . 

Therefore, 

3  lnp  8  h  (0  /Y)  M  9  lnp  y/9  (Y/e) 

90  0  -  0ML(Y)  90  ©  -  0  ML(Y) 

This  is  refered  to  as  the  Maximum  Log  Likelihood  Estimation 

of  6 

9 lip y/6  (Y/0)  ■=  0 

90  0  -  e^CY) 

This  estimation  of  0  will  now  lend  itself  to  a  form  of  tracking 
loop  as  demonstrated  by  Verterbi  (Reference  11). 

The  probability  density  of  the  noise,  (Ni(t),  i  =  E,  &)  is 
known,  (independent,  white,  zero  mean,  gaussian  with  identical 
independent  variances  V),  and  therefore,  the  joint  probability 
density  pf  (XL  (t),  XA(t)/a,e)  (f  (3fe  (t) ,  (t)/A,0)  can  be  derived. 

The  Gram  Schmidt  or  logonalization  procedure  (Reference  8) 
was  used  for  this  derivation. 

Gram  Schmidt  OrthoRonalization  Procedure 

As  Ni(t),  i  m  E,A),  are  independent  white  gaussian  noise 
with  independent  identical  variances  then  the  noise  can  be  represented 
(Reference  8)  as; 

n(t)  -  2  NiH(t) 
i  -  1 

Where  <t>i(t)  is  any  complete  orthonormal  set  over  the  interval 
0,T  and  the  Ni  are  independent  random  variables  with  variance  V 
as  Ni(t),  (i  *  E,A)  are  white  and  independent  with  independent 
variances 
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Also  expanding  x(t)  on  the  same  basis  set 


X(t)  -  IXi4>i(t)  for  0  1  t  <_  T 
i-1 

and  as  XT  (t)  -  AgI(9)S(t)  +  NS  (t) 

XA(t)  -  AgA(8)S(t)  +  NA(t)  then; 

XT  i  *  V  T/* 

s  XI  (t)$i(t)dt  and  NSi  «  ^  NE(t)^i(t)dt 


XAi 


y 

'A 


XA  (t)<(> i(t)dt  and  NAi  =  NA  (t)4  i(t)dt 


Let  <(> x<t)  *  S(t) 

fJT 


E  -  B2T 


B  *  1 


Kti 

T 


XI 


1m  ^  AS(t)gI (9)S(t)dt  +  £  NI(t)S(t)dt 


/E 


/F 


Agl  (9  ) 


y  s2(t)dt  +  NI. 
EX)  L 


=  Agl  (8 )  E  +  NI, 


Similarly  for  XA, 


XA 


T  T 

x  -  ^AS(t)gA(8)S(t)dt  +  y  NA(t)  S(t)dt 


/e 


/e 


(8) 


XAX  =  AgA(0)/I  +  NAX  (9) 

Statistics  of  XI  ^  and  XA^ 

XI ^  and  XA^  are  both  the  sums  of  deterministic  signals  and 
gaussian  random  variables.  Consequently,  both  XI  ^  and  XA^^  are 
gaussian. 


S(t)  dt  =  0 

ft 


Var 


Similarly 


Var 


Var 


[NA(t)] 


NZ(t)NZ(j)  S(t)  S(/)  dt 

ft  ft 


<*<r 


] 


=  I2  y  s2(t)dt 

0 

No 

2 


TnA.I  =  e[  Y  NA(t)N&(jr)  S(t)  S(/)  dt  d/1 

L  °  «  ft  ft  J 


No 

2 

No 

2 


E  [xZ1/0,Aj 
E  [xA1/0,a] 

[xZ1/0,a] 


y  S2(t)dt 


E  [AgE(e)//E] 
Agz(e)  ft 
E  [AgA(e)  /Ej 
AgA  (8 )  ft 


+  E 


+  E 


[KrJ 

H 


=  E 
=  E 


(Ags(e)  ft  +  NZ1)]  -  (Agz(e)  ft)2 
(AgE(e)  ft)2  +  2AgZ(9)  /E<NZ1)2 
(Agz(e)  ft)2 

-  (Agz(e) /e)2  +  —  -  (Ague)  ft)2 


(10) 


(ID 


r 


Similarly  Var  XAi/0,A  - 

XX  -  AgZ(0)/T  +  NX 
XA  =  AgA(©)/¥  +  NA 

(At  this  stage,  analysis  will  refer  to  first  pulse  only.) 

Where  the  NZ  and  NA  are  independent  and  identically  distributed 

No 

gaussian  random  variables  with  mean  zero  and  variance  y- .  The  joint 
density  function  (f(XZ,  XA)/0,A)  is  given  by 

p(f(XZ,  XA)/0,A)  -  *  (Xi/0,A)  i  =  Z,A 

r— 

l 


No 
2  IT 


1  exp 


(XZ-AgZ (0)  r7  E)  2  (XA-AgA (9)  A E )2 

~  -No 

27- 


2n° 

2 


=  K  exp 
Maximum  Likelihood  Estimator 


(XZ-AgZ (8)  /E)2  +  (XA-AgA (9)  /E)2 

No 


p  (f (XZ ,  XA)/0,A)  can  be  maximized  by  setting  the  partial 

derivative  ofp(f(XI,  XA)/0,A)  with  respect  to  0  to  zero. 

The  partial  derivative  of  (f(XZ,  XA)/0,A)  will  be  aided  if 

first  the  Inp  (f(XZ,  XA)/G,A)  is  found. 

Inp (f (XZ ,  XA)/0,A)  =  InK 

(XZ-AgZ (9)  /E)2  +  (XA-AgA (6)  /I)2 

No 

InK  -  XZ2-2XZArZ(6)  /!  +  (ArZ(9)/e)2 

No 

XA2-2XAAgZ(9)  /E  +  (AgA (9 )  /l)2 
No 

lnK-l  XZ2+XA2-2A/E(gZ(0)XZ+gA(0)XA)  +  EA2(gZ2(0)+gA2(0) )  (12) 

NO 


||lnp(f(XZ,  XA)/A ,©) 


■No[_2A  E(XZgZ(0)  +  XAgA(9))  + 


2EA  (gZ  (0)gZ  (0)  +  gA  (0)gA  (9)), 


\ 
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-  [-  (XEgE(6)+XAgA(0))  +  EA(gI(0)gI(e)  + 

gA(0)  gA(0))] 

Set 

fg  lnp  (f  (XT ,  XA)/A,8)  -  0 
Therefore, 

fj  lnp  (f(XE,  XA)/A,0)  -  -  (XEgE(0)  +  XAgA(0))  + 

EA(gE(e)gE (9)+gA(0)gA(0)) 

=  0 


Therefore; 

XEgE(0)  +  XAgA(e)  =  EA(gE (0) gE (0 )+gA (0)gA (0) )  (13) 

and 

fjf  lnp  (f  (XE ,  XA) /A,0)  -  -  ~  [2 /E(gE(0)XE+  gA(0)XA)  + 

2EA(gE2(0)  +  gA2(0))J 

Settlng  |j  lnp  (f (XE ,XA)/A,0)  *=  0 
Therefore; 

2  E(gE(0)XE+  gA(e)XA)  =  2  EA  (gE2(0)  +  gA2  (0)) 


(gE  (8)XE+gA  (8)XA) 
E(gE2(0)  +  gA  2 ( 0 ) ) 


(14) 


Substituting  equation  14  into  equation  13; 
XZgE (9)  +  XAgA(e)  - 


(gE  (8)XE+gA  (6)XA)(gE(0)gE  (0)+gA(0)gA(0)) 
(gE2(0)  +  gA2(0)) 


|xEgE(0)+XAgA(0j  •  [gE2(0)+gA2(0)J  -  [gE (0)XE+gA (0)Xa| 

gA(0)gA(e) 


•  gE (9)gE (0) 
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Expanding; 

[gE2(e)xEgE(e)  +  ga(0)xAgE(0)gr(e)  +  gi(0)xigA(0)gA(0)  + 

g2A(0)gA(0)XA J  - 

[gE2(0)gE(0)XE  +  gA(0)XAgI2  (0)  +  g2A(0)gI(0)XI+^A(0)g2A(0)XA] 
Therefore,  after  cancellation; 

gA(0)gE(0)gE(0)XA  +  gA(0)gE(0)gA(0)XI  -  gA (0)g2I (0)XA  + 

g2A(0)gZ(0)XI 


Collecting  like  terms; 

CgA(0)gE(0)gA(0)  -  gA2(0)gZ(0))XZ  -  -  (gI(0)g£(0)gA(0)  - 

gE2(0)gA(0))XA 

gA(e)(gA(e)gz(0)-gA(0)gZ(0))xz  -  -  gZ(e)(gz(e)gA(e)  - 


gZ(0)  gA (0) )XA 


gA(0)XZ  *  -  gI(0)XA 

gA ( 0 )  XI  gl (0;XA  ;  (^(Qjg A(0)-gE(0)iA(0)) 


Therefore; 


gA(0)XI  -  gI(0)XA 

and  for  the  maximum  likelihood  estimate  of  0 

gA(0)XI1  -  gK^XAj^  *=  0  (15) 

From  equations  8  and  9 

XI x  -  Agl  (0)  /E  +  NZ 

XAX  «=  AgA(0)  ft  +  NA 

and 

XI  -  y  S(t)  XI (t)  dt 

°  ft 


XA,  - 


y  sitx  xa(o  dt 

i  ft 
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This  is  equivalent  to  performing  the  estimation  process  after 


the  incoming  signals  have  been  coherently  demodulated. 

Error  Voltage  e  produced  by  the  maximum  likelihood  equation 

A 

Equation  15  states  that  when  0^  is  equal  to  0  (angle  off 

boresight  of  the  target)  then; 

gA(0)XZ,  -  gE (0)XA-  is  equal  to  zero.  When  0„,  does  not 
1  1  ML 

equal  0  then  equation  15  will  equal  some  number  other  than  zero. 

This  number  is  the  error  voltage  and  through  the  discriminator 
(e  to  0  mapping)  a  0  estimate  of  the  angle  off  boresight  can  be 
made.  Setting  equation  15  to  e 

gA(0)XS1  -  gE (0)XA  =  e  (15a) 

where  e  =  0  when  0=0 
e  ^  0  when  0  f  0 


Statistics  of  e 

E  [e]  =  gA(0)  E  [xrf  -  gE  (0)  E  j>4) 

from  equations  10  and  11 

E  [ej  =  gA (% ) AgE  ( 0 )  /E  -  gE  (§)AgA  (0)  /E 
E  [e]  =  A/E(gA(§)gE(0)  -  gE(@)gA(6)) 


If  6  =  0  then 
E 


-2 


[«]  =  0  <16> 
Var[e]-E  [(gA(0)XE-gE(0)XA)2]  -  A2E(gA (0)gE (0)-gE (0)gA (0) ) 2 

-  E  [gA2(0)XE2  -  2gA(0)gE(0)XEXA+  g2E(0)XA2]  -  e 
-[EA2gA2(0)gE2(0)  +  gA2(0)  +|2gZ2(0)  - 

2gA (0)gE (0)E  [xEXa]  -  e~2  +  EA2g2E (0)gA2(0)| 

(AgE (0)  /E  +  NE) (AgA (0)  /I  +  NA) 


[xZXa) 


-  E 

-  E 


'>] 


A  EgE  (0)gA  (0)  +  AgE  (0)  /E  NA  +  AgA(0)/-rE  NE  + 


NENA 


] 


! 


-  A2EgE(0)gA(0)  +  E  [nZNa] 

after  cancellation  of  zero  terms;  AgE(0)  E  NA,  AgA(0)  E  NT  and 
E  [NINA] 

Var  [e  ]  *  EA2gZ2(0)gA2(0)  +  |^[gA2(0)  +  gE2(0)]  - 

2EA2gI(0)gA(0)gA(0)gZ(0)  - 

E"2  +  EA2gA2(0)gZ2(0) 


EA2  JgE2(e)gA2(0)  +  gE2(0)gA2(6)J  -  2  EA2gE  (9)gA(e)gA (0)gZ  (0) 

No  f 

rlg 


+  ^|gA2(6)  +  gE2(9)j  - 


e“2  =  A2E(gA (0)gE (0)-gE (0)gA(0))2 

=  A2E(gA2(0)gZ2(0)  -  2gA(9) gE (0)gE(0)gA(9)  +  gE2(0)gA2(0)) 

e"2  *  -A2E(gA2(e)gE2(0)  =  gE2 (9)gA2 (0) )  +  2A2EgA (6)gl (0)gE 

(6)gA(0) 


Therefore  other  cancellation; 

N°  [gA2(0)  +  gE2(0)] 


Var  e  = 


2  /  •  s4"  vv/iJ  (17) 

The  density  function  of  e  is  the  sum  of  two  gaussian  functions 
that  are  multiplied  by  constants  gA(0)  and  gE(0),  and  is  consequent¬ 
ly  gaussian  with  the  following  statistics  wi,th  0  =  0; 


=  0 


e2  =  |2-  [gA2(0>  +  gZ2(0)] 

Statistics  of  A  estimate 


(16) 


(17) 


from  Equation  14 

A 

A 


(gZQ)XE  +  gA(e)XA) 
/E  (gE2(0)  +  gA2(9)) 
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r 


e  [X]  +  i 


E(gz2(0)  +  gA2(e>) 


[gr(e)E  XE  +  gA(0)E 


let 


from 


K  - 


E(gE2(0)  +  gA2(6)) 


E  [l/^  -  K  [gE(0)AgE(0)  /E  +  gA (0) AgA (8) 

-  K  A  /t  [gE(0)  gE(0)  +  gE(0)gA(0)] 

E  .  A_£SE,ie^E(92__±_jS.L(9i^.leJ 


2  A  2  a 

gE  (6)  +  gAZ(0) 


Under  the  condition  of  0  ■  0  ;  then 
E  [a/a]  -  A 

A 

Therefore,  if  0  =  0  ,  then  the  estimate  of  A  is  conditionally 
unbiased. 

Var  [a/aI  -  E 


[gZ(0)XZ  +  g A(0)XA) 

2  /A 

E  (gZ2(0)  +  gA2(0) 

_ 

-A 


2  *  2  A  2 
E(gZZ(0)  +  gAZ(0))^ 


E  [(gI(8)XZ  +  gA (0)XA) 2/a]  -a' 


2  A  2  a  2 

E(gE  (0)  +  gA*(0 ) ; 


(18) 


(19) 


(g2E(0)X2E  +  2gZ(0)gA(6)XEXA  + 


gA2(0)XA2)/A 


-A 


let; 


o  a  2  a  2 

E(gZ/(9)  +  gA  (0)) 


Var  [a/a]  -  KE  (g2Z (8) (A2g2Z (0) E  +  2AgZ(0)/?  NZ  +  (NE)2)  + 

2gZ(0)gA(0)(AgZ(0)  ffE  +  NZ)  (AgA  (0)  /~E  +  NA)  + 
g2A(0)(A2g2A(0)E  +  2AgA (8 )/E  NA+  (NA)2]  -A2 
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-  k  [A2Eg2z(e)g2z(e)  +  g2z(e)|S  +  2gE(e)gA(e)A2Egz(e)gA(e)  + 

A2Eg2A(6)g2A(0)  +  g2A6~]  -A2 

r  K  [a2e  g2z(0)g2r(e)  +  g2A(e)g26(e)  +  2gz(e)gA(0)gz(0)gA(0)>+ 
f2-  (g2Z(0)  +  gA2(0))]  -A2 


again  let  0  ■  0  ; 


[a/a]  -  - 


-5 - = - ~  A2E(gZ2(0)  +  gA2(6))2  + 


E(gZ  (0)  +  gA  (0))- 


(g2Z(0)  +  gA2  (0) 2 


after  cancellation; 


Var  [a/a]  -  Is  - 


(gZ2(0)  +  gA2(0)) 


To  determine  if  the  conditionally  unbiased  estimate  of  A 
is  an  efficient  estimate  then  the  conditional  variance  equation  20 
needs  to  be  equal  to  an  established  lower  bound,  (Cramer-Rao  Bound, 
Reference  8). 

Cramer-Rao  Bound  (Reference  8  ) 

Var  [A/A]  >  (e  [-—[lnp  (f(XZ,  XA)/A6)]  fj)"1  ' 

.  ra2  ]  N-i 

°r  -  (-E— y  lnp  (f(XZ,  XA)/A6)  ) 

'  L3A  J  ' 

from  equation  12; 

lnp  (f(XZ,  XA)/0,A)  -  InK  -  No  XI2  +  XA2  -  2A  E  (gZ(6)XZ  + 

gA(e)XA)  +  EA2(gz2(e)  +  gA2(e))  ( 

lnp(f(XZ,XA)/e,A)  -  ^  2  E(gZ(0)XZ+  gA(e)XA)  + 

2EA(gZ2(0)  +  gA2(0)) 
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pulse.  The  analysis  will  now  be  extended  to  two  pulses  as  in 
Figure  9 . 


Figure  9.  Two  Receive  Pulses 

xx(t)  -  Agz(e)s1(t)  +  Agi(e)s2(t)  +  Ns(t) 

XA(t)  -  AgA(9)S1(t)  +  AgA(8)S2(t)  +  NX(t) 

Using  Gram  Schmidt  Orthogonalization  procedure,  decompose  the  input 


signals  into  two  dimensional  random  variables. 


xx(t)  =  xi^s^t)  +  xs2s2*2(t> 

XA(t)  =  XAjS^t)  +  XA2S2^2(t) 


Where  XX  -  (XX^  XX  £) 
XA  -  (XA^  XA2) 
S  -  (Sl>  S2> 
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Given  that  0(t)  is  constant  during  pulse  and  during  the  pulse 
repetition  interval. 


Therefore, 


|iM.  „ 


Then  (9)  will  not  change  during  the  pulse  or  from  pulse  to  pulse 
so  g^(9)  (i*£,A)  will  remain  constant  from  pulse  to  pulse  for  Z 
channel . 

S  -  AgS(e)  y  S(t)4>.  (t)dt 
1  0  L 

-  Agi(e)  y  s(t)  ~  dt 
=  AgE(e)/E 

S2  =  AgZ (0)  y  S(t)^2(t)dt 

=  Agi(e)  T  s(t)  ~  dt 
0  E 

=  Agr(e)  /¥ 

Similarly  for  A  channel 

s1  =  AgA(e)  /¥ 

S2  =  AgA(0)/E 

Therefore, 

XE  -  (AgZ(0)/E  +  NE  ,  Agl  (0)  /E  +  NT,,) 

XA  -  (AgA (0)  /~E  +  NA  ,  AgA(0)/E  +  NA2) 


Now  to  write  (XE ,  XA/9.A) 

Since  NT  and  NA  are  independent 
‘{((XE,  XA/0,A)  -  £(XE/0,A)  (XA/0A) 
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and  since  NZ^  and  NI^  are  independent  and  NA^  and  NA2  are 
independent ; 


MP,  XA/Q.A)  -  KXI^G.A)  f>(XZ2/0,A)  KXAj/G.A)  l>(XA2/0,A) 
therefore, 

t>(XE,  XA/0.A)  -  K  exp 

(XAj-AgA(0)  E)2  +  (XA2-AgA(6) /i)2 


I'k  +  (XZ2-AgE(0)//E)2 


lnp  (Xz,  XA/0.A)  *  InK  -  ~ 

—  —  NO 


(Xlj-Agzte)  /E)2  +  (XZ2-AgE(e) /e)j 


+  (XA^AgACe)  /e)2  +  (XA2-AgA(0)/E): 


(23) 


A 

For  Estimate  of  A 

2(xz1  -  Agz(e)  e) (-)gz(e)  /e  + 
2(xz2  -  Agz(e)  E)(-)gz(e)  /e  + 


3  In 


~p  (XZ,  XA/O.A) 


1_ 

No 


2(XAX  -  AgA(0/  E) (-)gA (0)  /I  + 


2(XA2  -  AgA(0)  E)(-)gA(0)  /E 
setting  lnp  Qi?*  °  0 

0  -  -XZjgE^)/"!  +  AgZ2(0)E  -  XZ?gZ(0)  /I  +  AgZ2(0)E  - 


-  XA  gA(0)  /I  +  Ag2A(0)E  -  XA  gA(0)  /E  +  AgA2(0)E 

A  T~  ^  —i 

2AE(gs2(e)+gA2(0))*/E 


gZ  (0)  (XZx+XZ2)+gA  (0)  (XAx+XA2) 


Therefore; 

l  m  gZ(0)(XZ1+XZ2)  +  gA(0)(A1+XA2) 
2  E(gZ2(0)  +  gA2(0)) 


(24) 
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For  Estimate  of  6 

From  equation  23; 

2(XEa  -  AgZ(e) 

2(Xr2  -  AgE(e) 

2(XAX  -  AgA(e) 

2(XA2  -  AgA(e) 

setting  -|^lnp(X£,  XA/0, A)  =  0; 

o  -  -a/e  xz^zce)  +  a2e  gr(e)gEe- 

-a/e  xz2g5;(e)  +  a2e  gi(e)gze- 

-a/e  XA^aCe)  +  a2e  gA(e)gA(e)- 

-a/e  x&2gA(0)  +  a2e  gA(e)gA(e) 

o=  -gz(e)(xr1+x5:2)-gA(e)(XA1+XA2)  + 

2A  E(gE(0)gE(0)+gA(0)gA(0)) 
using  equation  24; 

2  /e  [gE(0)xE1+X2:2)  +«A(e)(XA1+XA2)]»[gE(e)gr(e)+gA(e)gA(e)J  - 

2/1  (gZ2(e)  +  gA  2  ( 0 ) 

gKeXXTj^+XSg)  +  gA(0)(XA1+XA2)  (25) 

Cross  multiplying; 

gE2(0)gE(0)(X2:i+)a:2)  +  gA(0)gE(0)gr(0)(XA1+XA2)  + 

L 

gE(0)gA(0)gA(0)(XI1+  XZ2)  +  gA2(0)gA(0)(XA1+XA2)  - 


|g  lnp  (XI,  XA/©,A)  -  -  ~ 


E)(-)AgE(0)  /E  + 
E)(-)AgE(0)  /E  + 
E)C-)AgA(e)  /E  + 
E)  (-) AgA(0)  /e! 
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WgzoxxZj+xzp  +  gi  (eigAOXXA^x^)  + 
_*2<0>s:(eXs:1+s:2>  +  gA2(e)gA(e)(XA1+XA2)j 
FgACejgKejIzceKXA^XAj)  +  gz(e)gA(e)gA(e)(xz 


gr2(e)gA(e)(XA1+XA2)  +  gA2(e)gz(e)  (xz 


pv 


1+  V 


1+X^j 


x&2)gz  (eXgA(e)s:  (e)-gz  (e)gA  (e))  + 


gz(e))  *  o 

. 


+(xz1+xj:2)gA(eXgz(0)gA(0)-gA  (0)gz(e)) 


(gA  (9  )  gZ  (0  )  -gZ  (0  )  gA  (9  )  )  (XA1+XA2)gZ(6)-(XE1+Xz2)gA(e)  -  0 


therefore. 


(X\1+XA2)gZ(0)-(XZ1+  xz2)g&(e)  =  0 


(XZ1+XZ2)gA(0)-(XA1+XA2)gZ(0)  =  0 


(K1+XZ2)g6(0)-Oai+XA2)gZ(0)  -  e 


(26a) 


Each  received  pulse  is  processed  individually  and  no  cross 
terms  exist.  The  analysis  of  the  single  pulse  consequently  holds 
for  more  than  one  pulse. 

The  effect  on  error  voltage  e  from  using  more  than  one  pulse 
Using  equation  26a; 


a'<S)-0Bi*V  .*<«> 


27 


e 


A  a  A  °2 

%xe  gA(0)  +  hxz9g&e - -  gze - -  gz(%)  - 

1  2  2 

^xr^ACe)  -  Jsx&1gE(§)  +  55Xi2gA(e)  -  %XA2gE(e)  -  e 
h  XEjgA^)  -  XAjgZ^)  +  h  XE2gA(e)  -  XA2gE(0)  =■  e 
from  equation  15a; 

he  x+  h  t2  m  e  (27) 

This  states  that  the  error  from  one  pulse  can  be  added  to  the  error 
from  the  second  pulse  and  the  final  error  averaged  to  obtain  a  new 
error.  As  the  mapping  of  error  e  to  estimate  angle  0  is  a  one  to  one 

A 

mapping!  each  0  estimate  obtained  for  each  pulse  can  be  added  to  the 

, 

next  0  estimate  for  the  next  pulse  and  then  averaged. 

Implementation  of  the  error  e  to  0  estimate  discrimination 
Using  equation  15; 

gA (0)  XE  +  gl(e)  XA  =  0  (15) 

Assuming  no  noise,  NE  =  Na  =  0; 

gA(0)AgE(0)  E  +  gE(0)AgA(6)  E  =  e 

Using  equations  6  and  7; 


sin  koa  sin  (0-  -j) 
a  d 

koa  sin  (0  -  -^) 


sin  koa  sin  (0  +  -jj)  -  sin  koa  sin(o 


-1)1 

_ r_ 


koa  sin  (0  +  2 ) 


koa  sin  (0  - 

(28) 


As  this  equation  is  a  combination  of  sine  functions  and  cannot  be 
reduced  effectively,  the  digital  computer  will  be  used  to  produce  a 
series  of  0  estimates  for  various  errors  e  .  The  error  voltage 
resulting  from  the  estimator  will  then  be  compared  to  the  various 


errors  and  a  corresponding  0  estimate  chosen.  This  estimate  will  then 
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•Ojcsl 


be  averaged  with  preceding  estimates  to  account  for  multiple  pulse 
processing. 

A 

From  equation  28  EA  scales  the  set  of  e/0  curves.  As  A  is 
unknown,  (A  is  contained  in  the  incoming  signal) ,  then  the  estimate 

A 

of  A,  equation  14,  is  used  to  estimate  A  and  is  used  to  scale  the 

A 

e/0  curves. 

From  equation  14; 


A 

A  = 


A 


A 


(gS(0)XS+gA(0)XA) 

E(g22(0)+gA2C0)) 


A  A 

A  is  a  function  of  both  0  and  the  amplitude  of  the  incoming 
A  A 

signal.  Only  when  0=0  does  A  =  A,  therefore,  the  scaling  of  the 

A 

e/0  curves  will  change  for  every  pulse. 

The  tracking  Loop 

The  system  at  this  stage  has  estimated  an  angle  6  from  one  or 
more  pulses.  From  equation  27  and  the  discussion  following  equation 


A 

27,  the  0  estimate  from  successive  pulses  could  be  averaged. 

The  tracking  loop  can  be  completed  in  two  methods  as 
follows: 


a. 


A 

50  pulses  can  be  used  to  estimate  50  0^  estimates. 
These  estimates  can  then  be  averaged  to  produce  a 

A 

0  estimate  angle  off  boresight,  (Figure  10)  as 
follows : 


A 

0  = 


50 


50 


i  = 


i 

1 


A 

The  0  estimate  can  then  be  used  to  steer  the  antenna. 
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Figure  10.  Implementation  of  System  (a) 

b.  Over  50  pulses,  each  pulse  can  be  used  to  derive  a 
a  a 

A  8  estimate  and  this  A0  estimate  averaged  with  the 

A 

preceding  A0  estimates.  This  average  A0  estimate  is 
then  added  to  the  current  estimate  to  form  a  new 
estimate  which  is  fed  back  to  the  estimate  antenna 

A  A 

functions  g£(0)  and  gA(0).  The  next  estimate  is 
then  derived  using  these  updated  antenna  functions, 
(Figure  11) . 

A0  i  =  —  (i-l)A8i_1+A9i  ;  i  =  1,2, ...50 

0±  =  J  (i-l)A0i_1+A,01  41_1  ;  i*l,2,...50 

50  1  ,  %  a  a  where  '  =  0 

0  -  —  j  (i-l)A01_1+Aei 


This  system  is  equivalent  to  a  tracking  loop 

A 

(Reference  11)  which  updates  the  estimate  of  6 
estimate  based  on  the  error  e  produced  from  the 
preceding  estimate. 


Figure  11.  Implementation  of  System  (b) 
Implementation  of  Tracking  Loo 


Using  equation  14,  a  digital  computer  program  was  written 
to  implement  each  block  of  Figure  12. 

The  error  e  voltage  was  determined  by  implementing  equation 
11a  directly  using  the  digital  computer.  The  pulse  width  T  was  taken 
to  be  20ms  and  the  number  of  noise  samples  from  a  gaussian  noise 

A 

generator  was  100.  The  first  estimate  value  of  0  used  in  the  feedback 
or  estimate  antenna  functions  gE(0)  and  gA(8),  was  arbitrarily  taken 
to  be  ziero.  The  error  voltage  e  was  then  compared  to  a  set  of 
numbers  generated  by  the  discriminator.  Each  of  the  discriminator 


3: 


bores ight. 

A 

The  probability  density  function  of  0  is  unknown.  Also  the 

A  A 

expression  for  0  is  unknown  as  0  is  contained  in  the  expression  for  e 
which  is  the  sum  of  sine  functions.  The  lower  bound  for  the  variance 
of  0/0  can  be  determined  using  the  Cramer-Rao  Bound  (Reference  8). 


Var  [e/e]  >  (e  [|^  In  p  ( f (XE ,  XA)/A0)  2]) 
>  ('E[ 


r-ln  p  (f(XE,  Xa)/A0) 


Var  [6/eJ  >  (-e(-§q  -  ~  -2A/E  (XEg(0)  +  XAgA(0)) 


2EA2  gE(8) gE (0)  +  gA(0)gA(0) 


I)'1 


>  [-E  [-  ^  -  2A  E  (XEgE (8)  +  XAgA(0))  + 

2EA2  (gE(0)gE(8)  +  (gE (0) ) 2  +  gA (0)gA (0)  + 

(gA(0))2)J  ^  -1 

Using  E  [xi]  =  AgE(0)/E  and  E  [xa]  =  AgE(0)/¥ 

Var  [e/e]  >_  +  ^  2  A2  E  (-gE (0)gE (0)-gA(0)gA(0)  + 


A  «« 


A  •« 


gE (0)gE  (0)  +  (gE (0) )  +  gA ( 0 ) gA ( 0 )  +  (gA (0) )  (34) 


Equation  34  was  derived  for  the  lower  bound  of  the  variance  of  0 

A 

or  the  best  possible  case  of  the  variance  of  8.  In  this  case,  the 

A  * 

best  possible  estimate  is  when  0=0.  Using  6=0. 


Var 


2-~ -  (-gE(8)gE  (8)-gA(6)gA(6)  + 


gE(0)gE(0)  +  gE2(0)  +  gA  (0)*gA  (0)  +  gA 2  ( 0  )  ) 


-1 
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TABLE  I.  Computer  runs  conducted  with  a  constant  target  at 
1.5°  angle  off  boreslght 


1st  Pulse 

50th  Pulse 

No 

2 

Var  j_6/1.5°J 

Var  6/1.5° 

.01 

1.0576  X  10~3 

-3 

2.1152  X  10-5 

.09 

9.5186  X  10  J 

-2 

1.904  X  10 

-4 

.25 

2.644  X  10 

-1 

5.288  X  10 

-3 

1 

1.05  X  10 

-1 

2.115  X  10 

-2 

4 

4.23  X  10 

-1 

8.46  X  10 

-2 

9 

9.5186  X  10  A 

1.904  X  10 

-2 

25 

50th  i 

pulse; 

"  ,  <3 

2.644 

5.288  X  10 

No 

V  1 

0/1. 5*] 

r  * 

50th  pulse 

2 

50EA2  X  4.7276  X  102 

No 

V 

[e/i.5 

| 

Nth  pulse 

2 

N  EA2  X  4.7276  X  102 

From  equation  35  the  lower  bound  of  the  variance  of  0 
estimate  is  a  function  of  the  noise  power,  the  amplitude  of  the 
incoming  signal  and  the  change  in  the  two  antenna  functions. 

In  the  computer  model  runs,  the  noise  power  and  the  amplitude 
of  the  incoming  signal  are  held  constant  for  each  individual  run 
but  in  the  case  of  a  moving  target,  the  antenna  functions  will 
change  with  change  in  angle.  This  will  cause  the  variance  estimate 
to  change  also  with  each  new  angle. 


Moving  Target 

All  theory  investigated  so  far  has  assumed  a  constant  target. 
However,  a  tracking  radar  needs  also  to  track  a  moving  target.  For 
the  computer  runs  with  a  moving  target,  the  target  will  assume  to 
move  from  -1.9°  to  3°  in  50  pulses. 

In  the  constant  target  model  progressive  averaging  of  the  A6 
was  done  over  the  50  pulses.  As  the  assumed  moving  target  is  almost 
covering  the  whole  radar's  angular  limits  within  50  pulses,  pro¬ 
gressive  averaging  over  the  entire  50  pulses  cannot  be  done.  To 
determine  how  many  pulses  to  average  and  what  type  of  averaging  is 
to  be  done  is  a  control  problem  and  is  not  within  the  scope  of 
this  thesis.  Consequently,  no  averaging  will  be  done  and  each 
estimate  will  be  fed  back  to  the  estimate  antenna  functions  to 
be  used  to  correlate  with  the  new  pulse.  The  investigation  of  the 
moving  target  situation  will  be  limited  to  investigating  the  radar's 
ability  to  track  a  moving  target  with  no  averaging  of  A0. 

A 

Output  of  the  tracking  loop  8 

A 

The  estimate  0  of  the  tracking  loop  was  investigated  in 
three  methods  for  each  of  the  target  models. 

For  the  constant  target,  the  three  methods  were  as  follows: 

a.  For  a  given  noise  seed  starting  value  the  estimate  0 
mean,  variance  and  (mean/variance)  square  was  plotted 
for  a  run  of  50  pulses.  This  was  repeated  for 
signal  to  noise  values  of  40db,  20db,  14  and  6db. 

b.  A  criteria  based  on  variance  and  (mean/variance) 
square  was  chosen  by  experiemnt  to  determine 
whether  lock  or  no  lock  was  achieved  during  the 
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50  pulse  run. 

/S 

c.  The  probability  density  function  of  0  was  plotted 
using  a  histogram  approach  (Reference  7) . 

For  the  moving  target,  the  three  methods  were  as  follows: 

a.  For  a  given  noise  seed  starting  value,  the 

A 

estimate  9  was  plotted  for  a  run  of  50  pulses. 

b.  The  probability  density  function  of  '0  was  plotted 
using  a  histogram  approach. 

c.  A  criteria  based  on  mean  square  error  of  the 
estimate  to  the  non-noise  tracking  value  was 
chosen  to  determine  lock  or  no  lock  was  achieved 
during  the  50  pulse  run. 

The  density  functions  of  the  constant  target  and  the 
moving  target  can  be  derived  in  many  ways.  Two  such  methods 
are  the  Histogram  Approach  (Reference  7)  and  the  Gram  Charlier 
Series  Approximation  (Reference  5).  However,  only  the  Histogram 
Approach  approximation  will  be  used. 

Histogram  Approach  (Reference  7) 

The  limits  of  0  were  chosen  to  be  -3°  and  3°.  This  line 
was  then  divided  into  30  points  or  0.2°  increments.  The  occurrences 
of  0  estimate  in  each  increment  was  then  counted  with  the  resultant 
number  divided  by  the  number  of  runs  (50). 


Therefore;  P  (9) 


N 

*i8Sci 


0 


■  cell  center  point  -  -2.9,  -2.7;;;2.9 
N  -  50 

oc^  ■  0;  -o.i  >_ e  -  e  i  >_  o.i 
1;  -0.1  <_  0  -  ^  <_  0.1 


mean 


Variance 


Vaif§J=  3J?  p(6)  e  2  -(E  0  )2 

i*l  1 

Gram-Charlier  Series  (Reference  5) 

When  many  random  variables  are  summed,  the  control  limit 
theorem  of  statistics  states  that  the  probability  density  of  the 
sum  random  variable  approaches  the  gaussian  distribution.  The 
Gram-Charlier  series  is  an  orthogonal  polynominal  expansion  whose 
first  term  is  the  gaussian  density  function.  In  this  case,  the 
second  order  of  approximation  was  used.  Accuracy  of  the  density 
function  is  increased  if  higher  orders  are  used. 


P  (9) 


1  £ 
°  i=0 


(i) 


A 

6  - 


Where  when  evaluated 


A 


d6 


C 

C 

c 


1 

3 


1  first  order 


-  —  ^  second  order 
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(34) 


2  is  the  3rd  central  moment  of  p(f)  normalized  by  oJ 
3  *  ■^r  /  (9  -  9  )  p(0)  d§ 

a> 

if  Mn  denotes  the  n  the  central  moment  of  P  (0 ) 

00 

Mn  “  /  6n  p(6)  d© 

then  ,  *=  m3  ~  3  +  2  mi 

o3_ 

and  H3  (1^1)  -  (i^l)3-  3(1^1) 

-  (KM)  -  ^  A  <*-ii)\(L=i) 

27 T 

Therefore  after  substitution  in  equation  34 

PC)  -  -i-  e-T  (VI)2  + 

rr  \  u 


+  i  i-  e  -  ULLlj((^l)3.3(isl^ 

~  3  2rr  \  > 


Ill  DIGITAL  COMPUTER  PROGRAM 


Digital  Computer  Programing  was  utilized  using  FORTRAN  5. 

The  main  flow  Flow  Chart  is  in  Annex  B.  The  basic  tracking  loop  was 
implemented  using  subroutines  called  from  Program  Main.  Variations 
in  the  output  were  then  conducted  within  Program  Main.  Listings  for 
the  basic  program  are  in  Annex  C. 

Basic  Tracking  Loop 

Ten  subroutines  were  used  to  implement  the  basic  tracking  loop. 
The  subroutine  ERID  was  used  as  a  mainline  calling  subroutine  to  the 
other  subroutines. 

Subroutine  FNE 

This  subroutine  implemented  the  basic  tracking  loop  with  calls 

to  all  the  relevant  subroutines  except  noise.  The  purpose  is  to 
a 

create  a  e  to  0  curve.  This  is  only  one  curve  which  was  created 
and  the  scaling  of  A  estimate  was  not  done  within  this  program.  With 

A 

each  change  in  9  estimate  within  the  tracking  loop  A  estimate  changes. 

A 

Therefore,  the  basic  e  to  9  curve  need  only  be  created  once  with  scaling 

a 

by  the  current  A  each  time  a  c  to  9  curve  is  required.  The  e  to  0 
curve  was  created  by  holding  the  input  angle  at  0°  and  varying  the 
9  estimate  from  -2.8°  to  2.8°  in  0.01°  steps.  This  created  a  set  of 
e  voltages  which  will  later  be  compared  with  the  e  from  the  actual 
tracking  loop  to  determine  a  9  estimate  which  would  cause  the  system 
to  have  0°  off  boresight. 

Subroutine  ERID 

Subroutine  ERID  implemented  the  basic  tracking  loop  with  calls 
to  subroutines  XPAT  and  ERROR.  The  basic  loop  included  noise  as  well 
as  some  defined  input  signal.  The  resultant  error  voltage  e  was  scaled 


AO 


by  A  estimate  and  then  compared  to  the  e  to  0  curves.  The 
resultant  6  was  then  returned  to  Program  Main  as  the  new  6  estimate 
to  be  used  on  the  next  pulse  processing. 

Subroutine  XPAT 

XPAT  created  the  XI  and  XA  signals  and  then  derived  the 
voltage  which  was  handed  back  to  Subroutine  ERID.  XPAT  called 
subroutines  NOISE,  SDPAT  and  ESTAN.  The  input  to  this  subroutine 
is  EST  which  is  the  current  value  of  0  estimate. 

Subroutine  ERROR 

The  e  produced  by  XPAT  was  scaled  by  A  estimate  in  ERID  and 

A 

then  compared  to  the  set  of  e  to  8  (curves  created  in  FNE)  in  ERROR. 
The  resultant  0  estimate  going  back  to  ERID. 

Subroutine  NOISE 

An  IMSL  routine  GGNML  is  called  from  this  subroutine.  GGNML 
creates  a  gaussian  vector  with  a  given  seed  starting  point.  The 
output  vector  (R(I))  is  gaussian  where 
R(I)  -  (Y(I)S  +  M 

where  Y(I)  is  gaussian  zero  mean 
unit  variance 
S  is  chosen  variance 
M  is  chosen  mean 

Subroutines  SDPAT,  SINC.  ESTAN  and  ESTDEG 

These  subroutines  implement  the  antenna  functions  gl(8)  and 

gA(6). 

Subroutine  MEVAR 

This  subroutine  was  used  to  derive  the  mean  and  variance  of  all 
past  9  estimates  for  any  given  run  of  50  pulses. 
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IV  RESULTS 


The  results  of  the  computer  program  runs  are  divided  into  two 
sections;  a  constant  target  with  50  pulses  are  sampled  and  progressive 
averaging  of  the  incremental  error  estimate  is  done  for  each  pulse, 
and,  a  moving  target  with  constant  velocity  where  no  averaging  of  the 
incremental  error  is  done  over  a  run  of  50  pulses. 

Constant  Target 

A  plot  of  mean,  variance  and  (mean/variance)  square  was  plotted, 
(Figure  Dl,  D2,  D3,  D4) ,  for  signal  to  noise  ratios  of  20db,  14db, 
lOdb  and  6db.  These  plots  are  ’once  only’  runs  in  that  the  same  noise 
seed  starting  point  was  used  and  consequently  are  meant  to  be  typical 
examples  of  the  radar  model  tracking  in  various  noise  conditions. 

Each  of  these  plots  demonstrates  the  radar  model's  ability  (variance 
of  the  estimate)  to  estimate  the  angle  off  boresight  (1.5°  off  bore- 
sight),  of  the  target  for  the  various  signal  to  noise  conditions. 

The  mean  shows  the  progressive  mean  of  each  of  the  estimates  over  the 
50  pulses,  and  the  progressive  variance  of  the  estimates  over  the  50 
pulses,  and  the  progressive  variance  shows  the  relative  deviation  of 
each  of  the  estimates  for  pulses  1  to  50.  The  (mean/variance) 
square  plot  shows  the  relative  goodness  of  the  mean  and  variance 
combination.  Due  to  the  averaging  of  the  incremental  estimate  A9 
in  the  discriminator,  then  the  variance  should  be  expected  to  decrease 
for  each  pulse  and  asymptote  to  a  particular  value.  On  the  other  hand 
(mean/variance)  square  should  increase  for  each  pulse.  Both  of  these 
conditions  are  for  a  locked  system.  If  variance  increases  and 
(mean/variance)  square  decreases,  then  the  system  is  considered  to  be 
unlocked.  These  criteria  are  good  measures  for  lock  or  non-loch. 


conditions  when  the  system  is  working  well  above  threshold  or  well 
into  threshold.  Around  threshold  the  (mean/variance)  square  approaches 
1  for  the  locked  condition,  however  (mean/variance)  square  must  still 
increase  for  a  locked  system. 

From  Table  I  a  lower  bound  variance  is  possible  of  2.644  for 
a  signal  to  noise  ratio  of  6db.  However,  an  accuracy  specification 
for  the  system  needs  to  be  assumed.  Assuming  an  accuracy  of  +  0.5°, 
then  the  lower  bound  on  the  variance  needs  to  1.  This  equates  to  a 
signal  to  noise  ratio  of  lOdb. 

From  the  plots  of  mean,  variance  and  (mean/variance)  square,  the 
20  and  14db  runs  (Figure  Dl,  D2)  were  both  lock  situations  but  the 
lOdb  (Figure  D3)  and  the  6db  (Figure  D4)  resulted  in  non-lock.  The 
6db  run  broke  lock  on  the  fourth  pulse  and  did  not  demonstrate  any 
indications  of  regaining  lock.  The  lOdb  (Figure  D3)  run  did  not 
break  lock  until  the  15th  pulse.  Although  the  mean  did  track  the 
target,  the  increasing  variance  indicated  that  if  the  system  was 
allowed  to  continue  for  more  than  the  50  pulses,  the  mean  would  also 
break  away  from  the  tracking  value.  The  variance  tended  to  oscillate 
with  the  mean  of  the  variance  increasing.  Also,  the  value  of  (mean/ 
variance)  square  had  started  to  decrease  suddenly  by  the  35' h  pulse. 
Under  the  criteria  of  increasing  variance  and  decreasing  (mean/ 
variance)  square  indicating  a  break  lock  condition,  the  lOdb  run  is 
considered  to  have  broken  lock.  However,  within  the  50  pulses  the 
mean  tracked  the  target  and  the  variance  increased  slowly  (1.5  over 
50  pulses).  This  run  is  an  indication  of  threshold. 

The  plots  of  probability  of  0  versus  0  (Figure  El  through  E24) 
display  runs  of  50  different  noise  seed  starting  points  for  signal 


Y 
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to  noise  ratios  of  40db,  20db,  14db  and  6db.  Each  plot  demonstrates 
the  probability  of  6  versus  6  at  pulses  of  1,  10,  20,  30,  40  and  50 
for  a  particular  signal  to  noise  ratio.  All  of  these  plots  were 
derived  by  the  Histogram  Approach  (Reference  7)  using  a  cell  width 
of  (.4°). 

The  40db  plots  (Figure  El  to  E6)  show  the  relative  frequency 
(out  of  50  runs)  of  a  particular  0  estimate.  Figure  El  shows  that 
all  50  runs  chose  an  ©  estimate  within  1.0°  and  1.4°  with  a  center 
or  mean  of  1.3°  as  all  values  are  within  one  cell  width  the  variance 
is  equal  to  zero.  On  the  tenth  pulse  (Figure  E2) ,  the  mean  has  now 
moved  to  1.56°  and  the  variance  is  8.4  X  10  This  indicates  an 
increase  in  variance  but  the  20,  30,  40  and  50th  pulses  (Figure  E3  to 
E6)  show  a  decrease  in  variance  until  the  50th  pulse  (Figure  E6)  shows 
a  mean  of  1.5°  and  a  variance  of  zero  (all  ©  estimates  of  the  individ¬ 
ual  runs  within  one  cell) .  This  agrees  with  the  computed  variance 
estimate  for  the  50th  pulse  of  2.1152  X  10  This  is  most  certainly 
within  the  .4°  cell  width. 

S' 

For  the  20db  plots  (Figure  E7  to  E12)  tracking  variance  close 


to  the  computed  variance  estimate  was  obtained. 
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TABLE  II. 


20db  Comparison  of  Computed  and  Actual  Var 


Pulse 

Mean 

Var  (V  o]  actual 

Var  |V  q'  Computed 

1 

1.308 

.0911 

0.1058 

10' 

1.36 

.0515 

0.0158 

20 

1.568 

.03937 

0.00529 

30 

1.668 

.0742 

0.00353 

40 

1.604 

.0516 

3.77  X  10“3 

50 

1.52 

.0612 

2.115  X  10"3 

The  pulse  1  plot  (Figure  E7)  indicated  a  variance  very  close 
to  the  computed  value  (difference  of  .0147).  As  the  number  of  pulses 
increases  (Figure  E12)  the  actual  variance  does  not  decrease  as 
quickly  as  the  computed  lower  bound  of  the  variance  estimate.  How¬ 
ever,  this  difference  can  be  partially  accounted  for  in  that  the 
computed  variance  estimate  is  a  discrete  lower  bound  estimate  where  as 
the  actual  estimate  is  based  on  a  cell  size  of  0.4°.  The  final  mean 
1.52  (Figure  E12)  and  variance  of  0.0612  gives  a  final  6  estimate  of 
1.27°  or  1.77°  which  is  well  within  the  tolerance  of  1.5°  +  0.5°. 

Again,  the  initial  variance  of  the  14db  plot  was  very  close 
to  the  computed  variance  estimate  of  Table  I. 

TABLE  III.  I4db  Comparison  of  Computed  and  Actual  Var  [*9/ s] 


Pulse 

Mean 

Var  [s/s]  Actual 

Var  [0/0]  Computed 

1 

1.432 

.4554 

.423 

10 

1.524 

.238 

.0423 

20 

1.668 

.181 

.0212 

30 

1.584 

.344 

.0141 

40 

1.56 

.216 

.0106 

50 


1.556 


26 


.0085 


r 


generated  numbers  corresponded  to  an  angle  off  boreslght.  Once  e 
was  found  to  be  within  the  limits  (equating  to  0  +  0.01°)  of  a 
discriminator  generated  number,  then  a  corresponding  angle  off 
boreslght  was  chosen.  This  angle  was  then  added  to  the  original 
estimate  of  zero  and  this  new  estimate  was  then  used  as  (6)  in  the 
estimate  antenna  functions  gE(0)  and  gA(^).  The  process  was  now 
repeated,  but  with  the  resulting  angle  off  boreslght  from  the 
discriminator  being  averaged  with  the  previous  angles  from  the 
discriminator  before  being  added  to  the  previous  or  existing  value 
of  0.  This  new  0  was  then  fed  back  and  used  to  update  the  estimate 
0  used  in  the  antenna  functions  g£(l)  and  gA(0). 

The  digital  computer  implementation  will  now  process  a  number 
of  returns  from  a  target  and  produce  an  estimate  value  of  0,  (averaged 
over  the  number  of  pulses  processed),  for  the  angle  the  target  is  off 
boreslght.  This  angle  of  0  can  now  be  used  to  update  the  angle  0 
in  the  antenna  functions  gl(0)  and  gA(0).  In  a  real  system,  a 
servo  motor  (as  opposed  to  phased  array)  is  usually  used  to  steer 
the  antenna.  The  effect  of  the  servo  is  to  average  all  incoming  esti¬ 
mates  with  a  turning  to  the  resultant  mean. 

Statistics  of  9 

The  problem  is  now  to  determine  how  many  pulses  to  process, 

A 

(mean  of  0  =  0)  and  to  what  degree  of  certainty  (variance  of  0),  that 

/V 

0*0. 

A 

0  is  derived  via  the  discriminator  from  e.  Under  noise 

A 

conditions,  the  mean  of  e  is  equal  to  zero  only  when  0  *  0.  As  e 

A  A 

and  6  is  a  one  to  one  mapping,  then  the  mean  of  0  in  the  limit  as  the 
number  of  pulses  processed  becomes  large  (infinity),  is  the  angle  off 
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As  the  number  of  pulses  increase,  the  variance  decreases  from 
pulse  1  to  pulse  10  but  then  oscillates  around  an  average  value  of  .28. 
This  is  a  very  good  indication  of  threshold  where  the  variance  is  not 
decreasing  significantly  but  also  is  not  increasing.  Plots  for 
pulses  1,  10,  20,  30,  40,  and  50  (Figure  E13  to  E18)  show  the  mean 
unchanged.  The  final  mean  1.556  and  variance  of  .26  produce  an 
estimate  8  ,  1.556  +  .509  which  is  just  outside  the  assumed  limit  of 
1.5  +  .5°  (2.06°  and  1.04°). 

As  can  be  expected,  the  6db  plot  (Figure  E19  to  E24)  demonstrated 
a  system  operating  beyond  threshold. 

TABLE  IV.  6db  Comparison  of  Computed  and  Actual  Var  f© /©I 


Pulse 

Mean 

Var  [e /f]  Actual 

Var  Computed 

1 

1.46 

1.44 

2.644 

10 

1.396 

1.25 

.2644 

20 

1.656 

1.71 

.1322 

30 

1.564 

1.26 

.088 

40 

1.504 

1.64 

.066 

50 

1.632 

.97 

.05288 

The  system  still  demonstrated  a  mean  tracking  of  the  target. 
From  Figures  E19  to  E24,  the  system  displayed  a  correlation  at  3°. 
This  was  caused  by  a  routine  within  the  model  which  limits  the  range 
of  the  system  to  +3°.  As  the  variance  of  this  system  was  large 
(average  of  1.38)  each  individual  pulse  numbers  standard  deviation 
took  the  estimate  to  the  outer  cell.  Consequently,  all  estimates 
which  exceeded  this  cell  were  lumped  into  its  histogram. 


To  determine  a  probability  of  break  lock,  an  assumption  needs 
to  be  made  as  to  exactly  what  constitutes  a  break  lock.  As  has  been 


shown,  the  system  can  still  track  the  target  but  have  out  of  limit 
estimates.  This  will  only  occur  when  the  out  of  limit  estimates  do 
not  exceed  the  discrimination  limits  of  +2.8°.  If  the  system  tracks 
a  target  but  the  estimate  is  out  of  limit,  then  the  system  is  not 
operating  as  specified,  and,  consequently  will  be  considered  to  have 
broken  lock. 

From  Tchebycheff  inequality  (Reference  9); 
p  (e  -  n)>  ka  £ 

assuming  as  before  k  =  0.5°  n  =  E 

o^  =  Var  6/0 

Table  V  gives  a  comparison  of  the  probability  that  the  system 
will  be  within  tolerance  (+0.5°)  at  the  50th  pulse.  Table  VI  gives 
a  comparison  of  the  probability  of  break  lock  of  the  system  at  the 
50th  pulse. 

2 

P(n  ~  e  <  6  <  h  +  e)>  ~ 

e 

TABLE  V.  Comparison  of  Probability  of  Lock 


s 

2 

n  -  .5  <  9  <  i 

N 

n 

a 

40db 

1.5 

2.1152  X  10"5 

1 

20db 

1.52 

0.0612 

.98 

14db 

1.556 

.26 

.73 

6db 

1.632 

.97 

0 
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TAVLE  VI.  Comparison  of  Probability  of  Break  Lock 


s 

N 

n 

Break  Lock  Probability 

40db  _ 

1.5 

0 

20db 

1.52 

.02 

14db 

1.556 

.27 

6db 

1.632 

1 

Moving  Target 

Plots  of  the  estimate  (Figure  FI  to  F4)  over  a  run  of  50 
pulses  for  signal  to  noise  conditions  of  no  noise,  40db,  14db  and 
6db.  These  plots  were  made  to  give  an  indication  of  the  tracking 
ability  of  the  system  for  a  moving  target  (-1.9°  to  3.0°  over  50 
pulses)  and  no  averaging  from  pulse  to  pulse.  As  was  seen  in  the 
constant  target  case,  the  system  should  track  the  no  noise  value 
of  the  estimate  on  the  average  but  in  the  moving  target  the  variance 
of  the  estimate  cannot  be  expected  to  decrease  with  the  number  of 
pulses  (no  averaging).  In  the  above  cases,  the  radar  tracked  the 
target  but  with  varying  degrees  of  accuracy . ^  This  accuracy  will  be 

A  /v 

determined  from  the  probability  of  0  estimate  versus  0  plots. 

The  40db  plot  (Figure  G1  to  G6)  shows  the  mean  changing 
from  -1.1°  to  2.9°.  The  maximum  variance  for  any  of  the  50  pulses  was 
0.009984  on  pulse  27,  This  equates  to  an  accuracy  of  0.1°  at  a  mean 
value  of  0.796°.  The  angle  for  the  non-noise  case  at  pulse  27  is  0.7. 
This  gives  an  error  at  most  of  0.196°  which  is  well  within  the  assumed 
accuracy  criteria  of  +0.5°. 

■The  14db  plots  (Figure  G7  to  G12)  show  the  mean  changing  from 
-1.196  to  2.752.  The  maximum  variance  for  any  of  the  50  pulses  was 
0.374  on  pulse  30.  This  equates  to  an  accuracy  of  +.61  at  a  mean 


value  of  1.028°.  The  angle  for  the  non-noise  case  at  pulse  30  is 
1°.  This  gives  an  error  at  most  of  .89°  which  is  outside  the 
assumed  accuracy  criteria  of  +0.5°. 

The  6db  plots  (Figure  G13  to  G18)  show  the  density  function  as 
a  relatively  uniform  function  (variance  greater  than  one) .  The  mean 
of  the  estimate  did  track  the  target  with  indications  of  -1.092°  to 
2.416°,  but  with  a  maximum  variance  of  2.243  on  pulse  17.  This 
equates  to  an  accuracy  of  +1.498°  at  a  mean  value  of  -.368°.  The 
angle  for  the  non-noise  case  at  pulse  17  is  0.3°.  This  gives  an 
error  at  most  of  -2.178°  which  is  well  outside  the  assumed  accuracy 
criteria  of  +0.5°.  Again,  correlation  occurred  at  the  3°  cell.  This 
causes  an  artificially  smaller  variance  at  these  angles. 

General 

In  all  cases  of  computing  accuracy,  the  control  problem  of  the 
system  has  been  overlooked.  This  problem  is  when  to  turn  the  antenna 
(how  many  pulses)  and  on  what  value  to  turn  the  antenna  (estimate  on 
some  form  of  mean  or  average) .  As  long  a  s  the  system  is  in  a  lock 
condition,  the  threshold  could  be  extended  by  using  an  average  or 
progressive  mean  in  the  constant  case.  In  the  tracking  case,  this 
problem  is  outside  the  scope  of  this  thesis. 

The  lock  no  lock  conditions  in  general  can  be  determined  from 
the  relative  scope  or  change  of  the  variance  or  (mean/variance)  square. 
If  the  variance  decreases  and  the  (mean/variance)  square  increases, 
then  the  system  is  locked.  For  the  moving  target  tracking  problem, 
no  relationship  can  be  determined  for  the  variance  as  the  number  of 
pulses  increases.  However,  the  lock,  non-lock  conditions  can  be 
determined  by  examining  the  mean  square  error  of  the  estimate  in 
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relation  to  the  non-noise  tracking  value.  If  this  exceeds  the 
specified  limits,  then  the  system  is  either  operating  under  too  low 
a  signal  to  noise  ratio  (in  our  case,  threshold  is  approximately 
20db)  or  the  system  has  broken  lock  and  is  no  longer  tracking. 

In  this  example,  the  specified  mean  square  error  limit  was  0.25 
(+0.52). 


51 


V  CONCLUSIONS  AND  RECOMMENDATIONS 
Theory  Conclusions 

The  estimation  system  is  to  be  implemented  to  model  the 
Amplitude  Comparison  Monopulse  Radar  is  as  follows; 


gA(0)  y X£(t)S(t)dt  -  g£(0)  y XA(t)S(t)dt  =  e 
0  0 
A 

Where  e  =  0  when  0=0 

The  amplitude  estimate  A  to  be  implemented  so  as  to  estimate 
the  amplitude  of  the  incoming  signal  is  as  follows; 


A  =  (gs(e)xz  +  ra(0)xa) 

E  (gZ2(0)  +  gA2(0)) 


A 

where ;  A 


Var  [a/a|  =  |£ 


gS2(0)  +  gA  2 ( 0 ) 


This  A  estimate  was  found  to  inversely  vary  the  e  to  0  curves 
of  the  discriminator. 

The  antenna  functions  to  be  implemented  were  found  to  be  as 
follows; 


gr(e)  = sln  [koa  sin  (9+  j>3 

koa  sin  (0+  ~) 

gA(e)  -  sin  [koa  sin  (0  +  |)] 
koa  sin  (9  +  Ij-) 


sin  (koa  sin  (0-  -|)j 

- j, - 

koa  sin  (0  -  -|) 

sin  j^koa  sin  (0  -  ^)] 
koa  sin  (0  -  ^ ) 


Where  koa  =  43.0003097 
|  =  1.2375° 
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The  A0  estimate  output  of  the  discriminator  can  be  progressively 
averaged  for  a  constant  target.  Averaging  other  than  progressive  aver¬ 
aging  over  the  entire  50  pulses  needs  to  be  done  for  a  moving  target. 

The  lower  bound  of  the  variance  estimate  for  the  1)  estimate  of 
a  constant  target  is  as  follows; 
r  a  "l  No 

Var  [e/ej  3  2EA2(gE2(0)  +  gA2(0) 

Results  Conclusions 

For  the  constant  target  with  progressive  averaging  of 
over  the  50  pulses,  the  threshold  was  found  to  be  determined  by  exam¬ 
ining  the  lock  condition  as  well  as  if  the  0  estimate  has  exceeded 
the  specified  accuracy  limits.  In  this  example,  threshold  was  found 
to  be  14db.  The  lock,  non-lock  criteria  for  this  system  is, 
a  ’lock'  is  determined  if  the  variance  decreases  and  the  (mean/ 
variance)  square  increases  over  the  50  pulses. 

For  the  moving  target  with  no  averaging  of  A§  over  the  50  pulses, 
the  threshold  was  determined  by  examining  each  individual  mean  square 
error  of  the  estimate  in  relation  to  the  non-noise  tracking  value. 

If  this  exceeds  the  specified  limits,  then  the  system  is  either 
operating  under  too  low  a  signal  to  noise  ratio  (threshold)  or  the 
system  has  broken  lock. 

Recommendations 

As  this  study  was  limited  in  scope  on  several  areas,  the 
following  areas  are  recommended  for  future  research; 

a.  Investigate  the  effects  of  varying  the  incoming  signal 
amplitude  A  on  the  threshold  and  break  lock  frequency. 

The  use  of  A  in  the  discrimination  curves  suggests 


that  a  saturation  level  of  A  will  cause  poor  tracking. 

b.  Reproduce  the  plots  of  probability  density  of  <5  versus 

A 

6  using  the  Gram-Charlier  Series  approximation.  This 
should  produce  a  more  accurate  determination  of 
threshold  than  the  Histogram  Approach  approximation. 

c.  Extend  the  system  to  two  dimensions  and  investigate 
the  effect  of  coupling  from  one  dimension  to  the 
other . 

d.  Investigate  the  control  problem  of  when  to  steer 
the  antenna.  The  type  and  amount  of  averaging 
needs  to  be  determined. 
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APPENDIX  A 


sin  142,993  sin  6  +  .92873  cos  6 


42.993  sin  0  +  .92873  cos  6 


sin  142.993  sin  6  -  .92873  cos  6 


42.993  sin  0  -  .92873  cos  0 

sin((42.993  sin  0)  cos  (.92873  cos0 


42.993  sin  0  +  .92873  cos  0 
sin((42.993  sin  0)  cos  (.92873  cos  0 


c 

os( (42.993  j 

1  sin  (.92873  cos  6)} 

42.993  sin0  - 

-  . 92873  cos  e 

The 

values  of  gl(0) 

and  gA(0)  depend  on 

the  value  of  0. 

In 

the  computer  runs  using  the  model  with  a  constant  target,  an 

angle 

off  boresight  of  1.5C 

> 

was 

used.  To  enable  the 

lower  bound  of 

Var 

e/0  = 

-  -O  . 

1.5  to  be 

determined,  the  four  terms  of  gl(0)  and 

gA(  0) 

need  to  be  evaluated 

at 

=  1.5°. 

1st  term 

-  8.553  X 

10~3 

+  9.012  -  5.506 

=  3.515 

2nd  term 

=  -3.055  X 

10~3 

-  15.12  -  3.5 

=  -18.636 

3rd  term 

-  8.916  X 

10~2 

+  56.302  -  599.00 

-  -542.6 

4th  term 

-  3.185  X 

10"  2 

+  157.626  +  382.14 

=  +539.79 

gi?(1.5°) 

=  321.2 

1st  term 

=  3.515 

2nd  term 

-  -18.636 

3rd  term 


+542.6 


Using  g£(1.5° 


(g£2(l.  5* 
From  equation 


Var  je 


4th  term  “  -539.79 

gA2(1.5°)  -  151.56 

and  gA(1.5°) 

)  +  gA2(1.5°))  *  472.76 

35 


> 


No _ 

2EA2  4.7276  X  102 
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APPENDIX  B 
Program  Flow  Diagram 


APPENDIX  C 


Program  Hating 


C 
C 
c 
c 
c 
c 
c 
c 
c 

PROGRAM  NAIM 

C  THIS  PROSRAH  SIMULATES  THE  SUFFICIENT  STASTII 

INTEGER  NR,  NA, MB,  NS 

REAL  ANGLE, THETA, GP AT ,GPATA,6PATB,GPAT1,6PAT2,PI,PHI 

1  ,DEG,9ETA,AHEAN ,BME  IN, VAR1,VAR2,SNC,STF , OCC (5C , 31)  ,SORT (33) 

2  , EST, VNULL, SBK1 , SUM  2,  EST 1, EST 2, A ,HBOX(561) ,MHU(5S) 

3  , ETA (561) ,POS »561), BOX  1561), DUHHISS) ,ALPHA(5  5I 

4  ,R( 1001 ,X1,X2,DME AN(53) ,BVAR(53) ,X( 33) , T  133) ,ORG 

5  , EHEAN, CVAR 

6  »SNR 

DOUBLE PRECISION  OSEED 

DATA  ANGLE,THCrA,GPAT,GPATA,GPATB,GPATl,GPArr2,ESr,»I 

1  , PHI, VNULL, SOHl, SUM  2,ESTl, EST  2, OEG, BETA , AMEAN,SMELN, VAR1 , WAP2 

2  ,SNC,STF,X1,K2,A, OPG 

C ALLPLOTS (0,,t,,9) 

CALL  PLOT<0.,-11.,*3) 

CALL  PLOT <2. ,5. 5,-3) 

NA«561 

N8*53 

NC*100 

DO  IBS  1*1,101 
R<I)*0.8 
100  CONTINUE 

00  106  1*1,50 
00  105  J*l* 30 
OCCtl, J)*0.0 
105  CONTINUE 
100  CONTINUE 

DO  110  1*1,53 
ALPHA ( I) *0*  0 
OUHM(I) *0*0 
DHEAN(I)*0.0 


a*********,***************,*********,*****,. 


PROGRAM  CONOEN 

THIS  PROGRAM  PLOTS  T HE  PR03ALITY 

OENSITT  FUNCTION  OF  THE  ANGLE 

OFF  BORESIGHT  FOR  A  CONSTANT  TARGET 
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BVAR(I)>0*0 

NHUIt»lil 

119  CONTINUE 

00  115  M|H 
SORT(X)»1.0 
Y(I»0.0 
YCD-0.0 
115  CONTINUE 

00  129  1-1,561 
ETA(X)«O.0 
pos(i>-o.o 
BOX (X)-0»0 
HBOXfDM.O 

120  CONTINUE  . 

0SEE0-126657.DI 

PI *3. 141592656 

CALL  FNE(PI,ANCLE,THE1A,GPAT ,GPATA ,GPATB, GPAT 1,G°AT2 
1  ,SUM1,SUH2«EST1,EST2,£TA,P0S,MB0X,NA) 

00  ISO  K-l, 50 
EST-0.9 
BETA-0.0 
STF-l.O 

SNR- (LO610 (1C1. 0/ST F**2. 01) *10.0 
DSEEO«OSEEO*FLDAT(X) 

DO  125  X—30,31,2 
SORT (11*32) /2) -FLOATCI)/1C.O 
125  CONTINUE 

00  120  X-1,30 
X ( X) -SORT (I ) 

SORT(I)«(SORT(I*1)-SORT(I))/2.0*SORT(I) 

120  CONTINUE 

CALL  ERIO (MR, ANGLE, THETA ,GPAT , 6PAT A ,GPATB, GPAT 1, *»  AT  2 

1  ,R,X1,X2,DSEED,EST, VNULL, PI, POS, BOX , PHI, ALPHA, ANE6N,  BMEAN 

2  , VAR1, 9AR2, OUMM, S NC ,STF, NA ,NB ,NC , BETA, DMEAN ,BVAR,6BDX,MHU> 
00  660  1-1,50 

IF  (OUHN(I)  *GE*  SORT (30))  THEN 
OCCCI, 30) *OCC( 1,30) *1*0 
GOTO  630 
ENDIF 

00  620  J-1,30 

IF  (OUNM(X)  *LE*  SORT(J) )  THEN 

OCC(I,J)*OCC(X,J)*1.9 

GOTO  631 

ENDIF 

620  CONTINUE 
630  CONTINUE 
660  CONTINUE 
130  CONTINUE 

00  600  L-1,50 
EHEAN-O.O 
CVAR-0.0 
00  671  H-1,30 

EHEAN-OCC (L ,M) /50.0*SORT(H)+EHEAN 
CVAR-OCC(L , M) /SO.Q* (SORT (N) **2. O+CVAR 
670  CONTINUE 

CVAR-CVAR-EMEAN**2. 0 


PRINT'</rT5,A,T2B,F5.2,r33,A//m3,A,T29,I2// 

,1  T15,A,F1B.2//T11,A,F19.2///T2,A/ 

2  T9,A/T9,A,T3*,F5. 2,T91, A/T9,A,T21,Z2//F> ’, 

3  'SIGNAL  TO  NOISE  RATIO  a »,SNR, *08 •, ’AT  PULSE’, 

A  (L) , 'MEAN  a »,EMEAN, 'VARIANCE  >’,CVAR, 

8  'FIG  PLOT  OF  THE  PROBALITV  DENSITY  FUNCTION', 

6  'OF  THETA  ESTIMATE  VERS  THETA’, 

7  'FOR  SIGNAL  TO  NOISE  RATIO  OF *,SNR, 'OB’, 

•  'AND  AT  PULSE', (L> 

X(31)a-3. 

T ( 31) at* 

X(12>a|./I, 

V(32>«l./9. 

00  955  J»l, 30 
V(J)aOCC(l,J)/5B.O 
955  CONTINUE 

CALL  AXXS<0.,I.,'THETA’,-5,6,,B,,X(31),X(32>) 

CALL  AXIS (3., I., 'PI THETA >»,♦•, 9., 90. ,Y (31) ,VC32)) 

CALL  PLOT (0,,9 ■ ,3) 

CALL  LZNE(X,V,3e, 1,1,0) 

CALL  PLOriO.,0.,3) 

CALL  PL0MS.,-L1.,~3> 

CALL  PLOT<B.8,l.5,-3> 

DO  960  la 1G, 50, IB 
DO  95B  J*l,  30 
Y<J)«OCC(I, JI/IO.O 
950  CONTINUE 

CALL  AXIS (O.,0., 'THETA', •5»6*,0.,X(31),X(32>> 

CALL  AXIS  C3«,0a,*P( THETA )',+B,9»,90,,Y(31) rY<32)l 
CALL  PLOMC.,1.,3) 

CALL  LINEtX.V, 30, 1,0,0) 

CALL  PLOT (B,,0* ,3) 

CALL  PLOTCB.,-11.,-3) 

CALL  PLOT (0*0,5 ,5, “3) 

96B  CONTINUE 

CALL  PLOTE(N) 

908  CONTINUE 
ENO 

SUBROUTINE  ERI D(NR,  ANGLE, THETA',’ GPAT,GPAT A, GPATB,SPAT1,  GPAT 2 

1  #R,X1,X2,DSEEC,EST,VNULL,PI,P0S,90X,PHI,ALPHA,AHEAV,BMEAN 

2  , VAR1, VAR2,DUHH,S  NC ,STF,NA,NB,NC,8ETA,DMBAN ,BVAR,NB0X ,MHU) 
REAL  POS(NA) ,B9X(NA), 0UHN( NB) , ALPHA (NB),MHU(NB) 

1  ,R(NC> ,DHEAN(NB) , BVAR(NB) ,MBOX(NA) 

CALL  XPAT (NR, A NGL E, THETA , GPAT, GPAT A ,GPATB, GPAT 1, SPAT 2 
1  ,R,X1,X2,DSEE0,EST, VNULL,PI,NC,STF, A) 

00  135  I«1,56D 

B0X(I)aNB0X(I)*((GPATl**2*GPAT  2**2) *2 ,0E*2) /( GPAf l'Xl»GPAT2'X2) 
135  CONTINUE 

CALL  ERROR (POS , BOX, PH 1,VNULL ,NA ) 

ANGLE-PHI 

CALL  ESTANi ANGLE, PI) 

ALPHA (2) a ANGLE 

EST-ANGLE 

OEG-EST 

CALL  ANOEG(DEG,PI) 
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0UNM(1)*0EG 
00  141  1*1,49 

-  cm  XPAT  ( NR, ARGL E» TH ETA  ,  GPAT  ,  GPAT A ,GP ATB  ,  GPAT 1, CPA  T  2 
1  ,R,Xi,X2,0$EED,EST,VNULL,PI,NC,STF, A) 

00  138  M*1,56D 

BOX(N)  -NB0X (N>  M  (GP  AT 1**2*GPAT2**2)  *2,0E-2)/(GPArt*Xl*6PAT2*X2) 
138  CONTZMUE 

CALL  ERROR(POS,BOX,PH),VNULL,NA) 

ANGLE*PHI 

CALL  EST1N (ANGLE, PI) 

ALPHA(I*2>*  ANCLE 
BETA*BETA*ALPHA  (1*2) 

EST*BETA/X*EST 

OEG-EST 

CALL  ANDEG(DES,PI) 

DUNN (1 +1) “DEG 
SNC*I+1 

CALL  MEVAR ( ALPHA, AMEAN*BME AN, VAR1,VAR2,SNC , DUNN, NB) 
0HEAN(X*1)-BHEAN 
BVAR(I»1>*VAR2 
HHU(X«1)*BHEAN»*2/VAR2*»2 
140  CONTI HUE 

RETURN 
ENO 

SUBROUTINE  FNE< PI, A NGLE, THETA, GPAT, GP ATA,GPATB,G>ATl,GPAT2 
1  ,SUN1,SUN2,EST1 ,E  ST2,ETA,P0S,MB0X ,NA) 

C  THIS  PROGRAM  ESTIMATES  T HE  PHASE  ERROR 

REAL  ETA(NA),POS(NA ) , NBOX(NA) 

ANGLE-0. 0 

CALL  EST AN (ANGLE, PI I 

CALL  SDPAT(ANGLE,THETA,GPAT,GPATA,GPATB,GPAT1,GPAT2) 

SUHlsGPATl 
SUM2-GPAT2 
00  220  X*.»28«,280 

ANGLE-FLOAT (I) (180.8  . 

P0S(X+281)*1 

CALL  ESTAN( ANGLE, PI) 

CALL  SOPAT (ANGLE, THETA, GPAT, GPATA, GPATB,GPATl ,GPAT2» 

EST1-GPAT1 

EST2-GPAT2 

ETA  (1*281  )*(SUN2*EST1-SUM1»EST2)*2.’0E>2 
220  CONTINUE 

00  230  1*1, 561 

NBOX (l)*(ETA(l*l) -ETA (I) )f  Z.O+ETAd ) 

230  CONTINUE 

00  240  1*1,561 
POSII)-P0S(X>/100,0 
240  CONTINUE 
RETURN 
ENO 

SUBROUTINE  ESTONIAN GL E,PI> 

C  THIS  PROGRAM  CONVERTS  DEG  TO  RAO 

ANGLE*AN6LE*PI/189 
RETURN 
ENO 

SUBROUTINE  SOPAT (ANGLE, THETA, GPAT , GPATA, GPA(TB,GPAT1,GPAT2) 

C  THIS  PROGRAM  CALCULATES  THE  SUM  AND  DIFFERENCE  AE  PATTERNS 
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THETA*AN6LE*.I2159BM»95 
CALL  PATT (THETA. GP AT) 

GPATA*GPAT 

THETA-ANGLE-. » 215986* 95 
CALL  PATT(THETA,GPAT> 

GPATB-GPAT 

GPATl*GPATA*GPATB 

GPAT2* GPATA-GPATB 

RETURN 

ENO 

SUBROUTINE  PATT (THETA, GPAT) 

IF  (THETA  .EG.  0)  THEN 

6PAT«1.I 

ELSE 

call  SXNC(THETA,GPAT) 

ENDXF 

RETURN 

ENO 

SUBROUTINE  SXNC (THETA. GPAT) 

THIS  PROGRAM  FINOS  GPAT 

GPAT*$XN(63<0I3B97*  SIN  (THETA)  )/  (43. B03B9M  SIN (THETA) ) 

RETURN 

ENO 

SUBROUTINE  XPA T (NR, ANGLE , T MET A , GPAT ,GPATA, GPATB,*PAT1, GPAT2 
1  ,R,Xl,X2, OSEEO, EST, VNULL ,PI,NC,STF, A) 

THIS  PROGRAM  CALCULATES  THE  XI  ,  X2  SIGNALS 
REAL  R(NC),STF 
A«1.B 
AN6LE«1.0 

CALL  ESTAN( ANGLE, PI) 

CALL  SOPAT (ANGLE, THETA, GPAT, GPATA,GPATB,GPAT1,GPAT2> 

CALL  NOISE (NR, R,OSEEO) 

Xl-B.t 

00  151  1*1, 1B0 

X1«((R(I)*STF*SPAT1)*A*2.BE-4)*X1 

CONTINUE 

CALL  NOISE (NR, R,DSEED) 

X2M.I 

00  161  1*1, 10 B 

X2*((R(I>  *STF*»PAT2 )• A*2  *BE-4) *X2 

CONTINUE 

ANGLE* EST 

CALL  SOPAT (ANGLE, TNET A,GPAT,6PATA,GPATB,6PATt,6PAT2l 

VNULL*X1*GPAT2»X2*GPAT1 

RETURN 

END 

SUBROUTINE  NOISE (NR ,R,OSEEDJ 

THIS  PROGRAM  GENERATES  A  NOISE  VECTOR 
REAL  A(ll(l 

DOUBLE  PRECISION  OSEEO 
NR* 100 

CALL  GCNNL  (OSEEO, nr, R) 

RETURN 

ENO 

SUBROUTINE  ERR  8  R  ( PO  S,  BOX  ,  PHI ,  VNULL  ,  NA ) 

REAL  POS(NA>,BBXir(A) 
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IF  (VNULL  .SC.  B0X(56O»  THEN 
PHI«P0S<»61> 

SOTO  CSS 
ENDIF 

00  CSS  X* It 566 

IF  CVNOLL  .LE.  BOXflll  THEN 

PHZ*POS(t> 

SOTO  Z6t 
ENOXF 

25  >  CONTINUE 
261  CONTINUE 
RETURN 
END 

SUBROUTINE  ANCEG<DEG,PX> 

C  THIS  PROGRAH  CONVERTS  RADIANS  TO  DEGREES 

0E6«0ES* IBS  /PI 
RETURN 
ENO 

SUBROUTINE  HEVARIAL PHA,AHEAN,8NEAN, VAR1.VAR2, SNC,  DUNH.NB) 
REAL  ALPHA (NB) .DUNN (NB) 

ANEAN«S.| 

8NEANM  .1 
VARl*t.B 
VAR2*a.s 
00  135  I* I. SNC 
ANEAN« ALPHA <X»2)«AMEAN 
VAR1*  ALPHA (I.21**2«VAR1 
BMEAN.OUNM<X>«BNEAN 
VAR2«0UNN(t)*»2«VAR2 
135  CONTINUE 

ANEAN«AMEAN/SNC 
VAR1«VAR1/SNC-ANEAN*»  2 

bhean-bnean/ss: 

VAR2«VAR2/SNC-BNEAN**2 

- -  RETURN . .  . . . . . ...  ,  .  ... - 

END 
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***»•»*••»•»*••»•#****»«♦*•♦***•***•**»**»••» 


c 

c  , 

C  *  PROGRAM  VEIHV  • 

C  •  • 

C  •  THIS  PIOGRAN  PLOTS  THE  ESTIMATE  » 

C  *  INGLE  OFF  BQRESIGHT  FOT  A  • 

C  •  MOVING  TARGET  * 

C  •  * 

c  •*••*»•****»***►•**•»***»****»*******»**•***» 

PROGRAM  MAIN 

C  THIS  PROGRAM  SIMULATES  THE  SUFFICIENT  STASH: 

INTEGER  NR»NA«MB«NS 

REAL  ANGLE»THETA» GPAT » GRATA fGPATBfGPATlt GP*T2» PI» PHI 

1  ,DEG,BETA,ANEAN,BME »Nf VAR1*VAR2,SNC»STF,VEL 

2  »EST#VNULL»SBMl,SUM2,ESTltEST2»A,MB0X(561),MHU(U> 

3  , ETA (5611  ,POS  (561)  (BOX  (561)  ,0UHH(53>  .  ALPHA  ( J  J> 

A  (R(10l)f  Xlt  X2f  ONE  AN  (S3)  (BVAR(53)  »X(53)  »  Y  (53) 

5  f  SNR 

DOUBLEPRECISION  DSEED 

DATA  AN6LE, THETA, GPAT, GRATA, 6PA TB, GPAT1,GPAT2, EST , PI 

1  f PHI* VVULLfSOMlf SUM2«EST 1*  EST  2>DEG(BETA , NME AN» BHE AN* V AR1, VAR2 

2  ,SNC,STE*X1,X2, A,  VEL 

3  •Qi0»0i0f0»0t0»Ct0tQf8fC>G|C«C(Ct(if9fttC/ 

C ALLPLOTS (0 •*  E • «9 ) 

CALL  PLOT(0.«-11>,*3) 

CALL  PLOT<E«(7*5(-S> 

NA-561 

NB*53 

NC-186 

oo  it  a  i«i«iat 

R<i>*o.o 
100  CONTINUE 

00  111  1*1(53 
ALPHA (I)*O.C 
0UMM(I)*».6 
DNEANTZXO.O 
8VAR(Z)*|.D 
MHU(I)«0.0 
X(II«0.0 
X(I>*0.0 
110  CONTINUE 

00  120  Z>1(561 
ETA (I)*0*0 
PDS(I)*0.0 
Box«)*o;o 

MBOX(I)*0»0 
120  CONTINUE 

DSEEO-12AA57.0I 
PI *3*161592651 

CALL  FNE(PI (ANGLE,! HE )A, GPAT, GRATA, GPATB,GPATl,G»AT2 
1  (SUM1(SUM2(EST1(EST2(ETA (POS( MBOX(NA) 

00  139  1*1(5 
EST*0.0 
BETA-0. • 

STF«I 

SNR* (L0G1I ( 101 «I/ST F**2*0) )*10*0 


59-6 


CALL  ERIOCNR, ANGLE, THETA, GPAT,GPATA,GPATB,GPAT1,;PAT2,VEL 

1  ,R,Xl,t2,0SEED,£ST,VNULL,PZ,P0S,B0X,PHZ , ALPHA, ANEAN,BHE AN 

2  , VAR1,VAR2,DUHM,SNC,STF,NA,NB,NC, BETA,DNEAN ,BVI R, (, NBOX, NHU) 
PRINT* C//T5 ,A,T28,F9.2,T33,A///T13, A,T35*I1,T44,I2,T54 ,I2S/ 

1  TU,A,Fia.2,Flt.2,F]9.2///T2,A/T9,A,Tir,F5.2,T22,l///>* 

3  , 'SIGNAL  TO  BOISE  RATIO  «*,SNR, «OB>, *AT  tt>  PULSES*, 1,2*, 56 , 

4  'ESTIMATE  -  * *,OUHH (1) ,OUMM (25> ,DUHH(50> * 

5  *FIC  PLOT  OF  ESTIMATE  FOR  SIGNAL  TO  NOISE* 

6  , 'RATIO  OF*, SNR, *06  AND  A  MOVING  TARGET.* 

130  CONTINUE 

ENO 

SUBROUTINE  ERIDIMR, ANGLE, THETA, GPAT,GPATA, GPATB,SPiri,GPAT2,VEL 

1  ,R,X1,K2,DSEED,EST , VNULL,PI,POS, BOX, PHI, ALPHA, AREA N, BHEAN 

2  ,VAR1,VAR2,D9MM,SNC,STF,NA,NB,NC,8ETA,DMEAN,BVAR,*,NB0X,NHU> 
REAL  POS(NA), BOX ( NA > , CUNH(NB) , ALPHA (NB) , HHU(NB) 

1  ,R(NC) ,ONEAN( NB) , BVAR(NB) ,X(NB> , HBOX(NA) 

VEL*1*B 

CALL  XPATINR, ANGLE, THETA ,GPAT,GPATA,GPATB,GPAT1,:PAT2 
1  , R, X1,X2,DSEED, EST, VNULL , PI, NC, STF, A , VEL > 

00  135  1*1,561 

BOX (I) *MBOX <1  )*  I  !GP AT  1**2*GPAT2**2)*2.|E“2)/I GPATI'X1*GPAT2*X2) 
135  CONTINUE 

CALL  ERROR! POS, BOX, PHI, VNULL, NA) 

ANGLE-PHI 

CALL  ESTAN!ANGLE,PI ) 

ALPHA  f 2)* ANGLE 

EST*ANGLE 

OEG*EST 

CALL  ANOEG(OEG,PI> 

0UMH(1)-DEG 
00  146  1*1,49 
VEL*FL0AT(I)*1.8 

CALL  XPATCMK, TOOLE, THETA, CPAT,GPATA,GPATB,GPATl,GPAT2 
1  ,R,Xl,X2,r%  it, ES T, VNULL, PI, NC, STF, A, VEL) 

DO  138  N*l,^l 

BOX (N)  ■MBOXfNi  *  ,  {(5PATl**2*GPAT2**2)*2.tE-2»/lGPAri»<l*GPAT2*X2) 
138  CONTINUE 

CALL  ERROR!POS,BCX, PHI, VNULL, MA) 

ANGLE*PHI 

CALL  EST1N(ANGLE,PI) 

ALPHA <I*2) -ANGLE 
EST*EST ♦ANGLE 
OEG-EST 

CALL  ANOEG (OEG, PI ) 

OUMMCIPI) *OEG 
SNC-I*1 

CALL  MEVARIALPHA,ANEAN,BMEAN,VAR1,VAR2,SNC^)UHN,V9) 

0MEANIIU1-BMEAN 

BVARd^l)  *VAR2 

HHU!I*1) »BMEAN**2/V  AR2**2 

X(I«1)*I»1 

CALL  LOO(T,OUHN,NB) 

X<51)«9. 

T!51)*-3. 

Xt52)*58./5. 

Y(92)*B.T4. 

CALL  AXISIB., I., 'NUMBER* ,-4,5. , B. , X (51) , X < 52) ) 
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CALL  'M.EAN  VAR  X  A.,  98.,  V  151 )  ,Y(52») 

CALL  PL0T(8.,I.»3) 

CALL  AXX8(5*,-2.,  *AEANSa/VAR»,-10,A.,9».tM5n,T(5en 
CALL  PL0T(a.,l.,3) 

CALL  PL0T(8.,-2.*-SI 
CALL  LXNE(K,Y»  58,1,8,81 
CALL  PLOT (S.,l«  ,31 
CALL  PLOT (10. 8, -11. i-3) 

-  CALL  PLOT (8 .*2. 5, -II  •  ■  -  r  r-  . 

IAS  CONTINUE 

CALL  PLOTE(N) 

RETURN 
END 

SUBROUTINE  FNECPI ,A NGLE.THET A, SPAT ,GP ATA, GPATB, GPATI, SPATE 
1  ,SUN1,SUM2,EST1,EST2,ETA,P0S,HB0X,NA>  I 
C  THXS  PROGRAM  ESTIMATES  THE  PHASE  ERROR 

REAL  ETA ( NA )f P9S( NA )fMBOX(NA) 

ANGLE*!. I 

CALL  ESTAN(ANGLE,PI> 

CALL  SOPAT(AN6LE»TH ETA »6PAT»6PATA, GPATB, GP ATI, GPAT2I 
SUM1*GPAT1 
SUH2*GPAT2 
00  220  I*-28G»268 
ANGLE*FLOAT (II  fils. 8 
P0S(I*2Bl)*: 

CALL  ESTAN(AN»LE»PXl  I 

CALL  SOPAT(ANGLE,THETA,GPAT,GPATA,GPATB,GP4ti,GPAT2> 

EST1-6PAT1 

EST2-GPAT2  < 

ETA(I*28l>»(SUN2*ESTl-SUMl*EST2l*2.0E-2 
22D  CONTINUE 

00  238  1*1,56# 

MBOX(I>«tETA(I*l)-ETA(I>l/2.8*ETA(Il 
238  CONTINUE 

DO  2A8  1*1,561 
POS(1)*POS(IISIQO.I 
2%0  CONTINUE 
RETURN 

END  1 

SUBROUTINE  ESTAN ( AN GLE.PIl  4 

C  THXS  PROGRAM  CONVERTS  DEG  TO  RAO 

ANGLE*AN6LE*PX  f 180 
RETURN 
END 

SUBROUTINE  SDPAT(ANGLE,THETA,GPAT,GPATA,GPAfTB,GPATl,GPAT2l 
C  THIS  PROGRAM  CALCULATES  THE  SUN  AMD|  01 FFERCN3E  AE  PATTERNS 

THETA«ANGLE*. 1 215984495 
CALL  PATT (THETA, GPATI 
GPATA*6PAT 

THETA* ANGLE-. 12159644 95 
CALL  PATT (THETA, GPATI 
GPATB*GPAT 
GPAT1«GPATA4€PATB 
GPAT2*GPATA-GPATB 
RETURN 


I 


■aft* 


—WO.I "  V 
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ENO 

SUBROUTINE  PATT (THETA, SPAT) 

'  IF  (THETA  .EQ.  0)  THEN 
GPAW.O 
ELSE 

CALL  SINC(THETA,GPAT) 

ENDXF 

RETURN 

END 

SUBROUTINE  SING (THETA, 6PAT) 

C  THIS  PROGRAM  FINOS  6PAT 

GPAT*SXN(I3. OK 3197* SIN (THETA) >S (*3.003097»SIN (THETII  > 

RETURN 
EH  D 

SUBROUTINE  XPAT (NR, ANGLE, THETA, SPAT ,GPATA,GPATB,6PITl, 6PAT2 
1  ,R,X1,K2, OSEED,EST, VNULL ,PI ,NC,STF, A ,VEL ) 

C  THIS  PROGRAM  CALCULATES  THE  Xt  ,  X2  SIGNALS 

REAL  R(NC) »STF 
A«i»0 

ANGLE*VEL/1B. 1-2.0 
CALL  ESTIN (ANGLE, PI) 

CALL  SOPAr(ANGLE,THETA,GPAT,6PATA,GPATB,GPATl,GPir2) 

CALL  NOISE(NR,t,OSEEO) 

Xt«0.0 

00  t5«  i«i,ne 

Xl*((R(I)*STF*GPATi)*A*2.0E-A)>Xt 
ISO  CONTINUE 

CALL  NOISE (NR, t,OSEEO) 

X2«0.l 

00  161  1*1,101  ( 
X2*((R(II*STF*CPAT2)*A*2.#E-A)*X2 
160  CONTINUE 
ANGLE«EST 

CALL  SOPAf (ANCLE, T(ET  I,GPAT,GPATA,GPATB,CP«T1,GPIT2I 
VNULL*X1»GPAT2-X2*GPAT1 
RETURN 
ENO 

SUBROUTINE  NOISE(NR,R,DSEED> 

C  THIS  PROGRAN  GENERATES  %  NOISE  VECTOR 

REAL  R(lll) 

DOUBLE  PRECISION  DSEEO 
NR* 110 

CALL  GSNNL  (DSEEO, NR, R) 

RETURN 
ENO 
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258 

268 


C 


135 


SUBROUTINE  ERROR  <PO  S, BOX , PHI , VNULL ,  HA > 

REAL  POS(NA) VB9X( NA > 

IF  (VNULL  .GE.  BOX(560»  THEN 
PHX«POS(f61> 

SOTO  268 
ENDZF 

DO  258  Z*lf 568 

IF  (VNULL  .LE.  BOX(D)  THEN 

PHX«POS(t) 

SOTO  268 

ENDIF 

CONTINUE 

CONTINUE 

RETURN 

ENO 

SUBROUTINE  AND E6 ( DE G, FI) 

THIS  PROSRAH  CONVERTS  RADIANS  TO  DECREES 
0EG*0E6*188/PX 
RETURN 
ENO 

SUBROUTINE  NEVAR (ALPHA ,ANEAN,BHEAN, VAR1, VAR2.SNC.  D JNH, NB) 
REAL  ALPHA (NB) , OUHN (NS) 

AHEAM8.8 

BHEANM.S 

VAR1-8.8 

VLROM 

00  135  I>1,SNC 

AHEAN-ALPHA (I»2)*AHEAN 

VAR1»  ALPHA (I»2)**2«VAR1 

BHEAN>OUNH(X)»BHEAN 

VAR2«DUHN(I)**2*VAR2 

CONTINUE 


AHEAN*AHEAN/SNC 
VARl-VARl/SNC-AHEAN**  2 

bhean>bhean/sn: 

VAR2*VAR2/SNC”BHEAH**2 

RETURN 

ENO 

SUBROUTINE  L00(  >» OUMH  iNB) 
REAL  V (NS) yOUMN (NB) 

00  278  I»l,50 

IF  (OUHH(I)  .LE.  -3 .8)  THEN 

OUHH(X)*-3.0 

ELSE 

OUHH(I)*OUHH(I) 

ENDIF 

IF  (OUHH(I)  .GE.  3.8)  THEN 

OUHH(I)«3.8 

ELSE 

OUHH(I)>OUNH(I) 

ENOIF 

r(I)OUHH(I) 

CONTINUE 


RETURN 
END  ' 
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»  » 

*  MO  SR*  N  CON  NV  * 

*  • 

•  THIS  PROGRAM  PLOTS  THE  MEAN  ,  VARIANCE,  » 

•  ANO  THE  (MEAN /VARIANCE) SQUARE  FOR  R  » 

•  CONSTANT  TARGET  * 

»  * 


*»***»»**,*»•*******»*•*»•••*****#**,»******» 

PROGRAM  MAIN 

C  THIS  PROGRAM  SIMULATES  THE  SUFFICIENT  STASH! 

INTEGER  NR,NA,NB,NC 

REAL  AN6l E, THETA, GP AT , GRATA, GPATB, GPAT1, GPAT2 , PI, PHI 

1  , DEG, BETA ,AMEAN,BMEAN>VAR1,VARZ,SNC,STF 

2  , ESTfVNUL L, SBHl,SUM2,ESTl,EST2,A,HBOXt561),HHUf5St 

3  ,ETA (B61) ,POS (561) ,B0X (561) , OUHMC53) VALPHAtl  J( 

A  ,R(  IBS)  , XI, X 2,  ONE  AN  153)  ,  BVARI 53)  ,X(53),Y  <53) 

5  ,SNR 

DOUBLEPRECISIOI  OSEEO 

DATA  *NGLE,THETA,GPAT,6PATA,GPATB,GPATl,GPAT2,ESr,PI 

1  ,PHI,  MULL,  S9M1,  SUM 2,EST1, EST2,DEG,BETA,NMEAN,BNEAN,  VAR1, VAR2 

2  ,SNC,STF,X1,(2, A 

3  /■,•,(, I, 1,1, 9,0, •,8,B,B»tyf,B,B,B,B,l,B,B,C,A,l,l,l/ 
CALLPLOTStO*,f.»9) 

CALL  PL0T(B.,-ll.,-9> 

CALL  PLOT (t«,7,S,-3) 

NA«561 

NB«5J 

NC»100 

oo  iso  x*i,ioe 

R(I)»OiO  . 

100  CONTINUE 

00  110  1*1,53 
ALPHA(I)-0.0 
0(/MNTI)«0.0 
OMEAN (I)*0,  0 
BVARf I)*9«9 

MHU(I)*0,0  » 

x(i)*o«o  * 

x(i)*e,o 

110  CONTINUE 

DO  120  1*1,561 
ETA(I)*9.0 
POS(I)*0»0 
BOX  (I)*0,0 
MBOX(I)*I.O 
120  CONTINUE 

DSEEO* 120657,01 
PI *3.161592656 

CALL  FMECPI, ANGLE, THE 1A,GPAT,GPATA, GPATB, GPAT1,G>AT2 
1  ,SUM1,SUM2,EST1,£ST  2,ETA ,POS,HBOX ,NA) 

X6JZY6B606SSH72 (3RX6 

EST-0.0 

BETAM.O 

STF*I 
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SNR*<ioGii<ioi.o/srr**2.i>>*ie,c 

CALL  ERID (NR, A  NGLE, THETA ,GPAT,GPAT  A,GPATB, GPAT1,»P AT2 

1  ,R,X1,X2,0SEED,EST, VNULL, PI, P0S.90X, PHI, ALPHA, AHEAN,BHEAN 

2  ,VARi,VAR2,DBHH,SNC,$TF, NA ,NB ,NC , BET A,OHBAN ,BVAR, X, NBOX, HHU ) 
CALL  LDO(Y,OH£AN,NB ) 

X(51) »8« 

r(ju«*j. 

X(52>-5B./5. 

Y(52)*6.F*. 

CALL  AXXS(A«,t,  , 'NUMBER' ,-6,5«,C»»X (51) ,X 152) ) 

CALL  AXXS(0m-2.»  'HEAR  VAR  X •  , *10,4. f90. , Y (E51 ) , V<52>  > 

call  PLor<a.fa««3> 

CALL  AXXS(5,,~2,,  ,MEANSQ/VAR*f“10»(»»»9D«*Y  tSl),Y(52>) 

CALL  PIOTU.,1.,3) 

CALL  PLOT (0«,-2»,-3) 

CALL  LlNE(X,V,S0,l,0,aT 
CALL  PLOMe.,1.,3) 

CALL  LOP(Y,8VAR,N8l 
CALL  LZNE(XfV,I8,l,5,*) 

CALL  PLOr<0.»Bt'3) 

CALL  LOQ(Y,MHU,NB) 

CALL  LXNE(X,Y, 58,1, 5,1) 

CALL  PLOT  (0  «  ,(*  •  ,3) 

CALL  PLOT(10.t,-U. ,-3) 

CALL  PLOT <8  •  ,7»  5,  -3  ) 

PRINT*(/'T6,A,r29,F5,2,T3*,A///Tl*,A,T36,Il,T*5,t2,r55,IZF/ 

1  Tl9,A,Fie.2,F10.2,F18.2//T15,A,F10*2,Fia*2,F10.2" 

2  T2,A,Fl0.2,FL0.2,Flt.2///T3,A,T3B,A/TlB,A,T36,F5.2,r*l,l///)' 

3  t 'SIGNAL  TO  NOISE  RATIO  *',SNR, 'OB', 'AT  €IJ  PU. SEi *, 2, 25, 5f , 

*  'MEAN  -  **,DNEAN(2) , 0MEAN(25),DNEAN(5B), ‘VARIANCE  X  »• 

5  ,BVAR(2),8VAR(25)  ,BVARC5B),MHEAN/VARIANOE>SQUARE  0  *• 

6  , MHU (2) »NHU(25),NHU(50),'FIG  PLOT  OF  HEAN,  VARIANCE  AND  <*, 

7  'MEAN /VARIANCE )SQU  ARE*, 'FOR  SIGNAL  TO  NOISE  RATIO  OF', SNR, 'OB' 
130  CONTINUE 

CALL  PLOTECN) 

ENO 

SUBROUTINE  ERIOINR, ANGLE, THETA, GPAT ,GPATA,GPATB,5PAT1, GPAT2 

1  , R,X1,X2,QSEED,EST, \NULL,PI,POS,BOX , PHI , ALPHA, AHEAN, BMEAN 

2  ,  VAR1,VAR2,0IIMM  , SNC  ,  STF ,  NA  ,NB  ,NC  »6ET  A,DNEAN  ,BVAR,X  ,  N80X,  HHU ) 
REAL  POS(NA),BOX(NA),CUHM(NB),ALPHA(NB),'NHU<NB) 

1  ,R(NC),DMEAN(NB) , BVAR(NB) ,X(NB) ,HBOX (NA) 

CALL  XPAT (NR, ANGLE, THETA,GPAT*GPATA,GPATB,GPAT1»»PAT2 
1  ,R,X1,X2, OSEEO,EST, VNU  LL ,PI,NC,STF,A) 

00  135  1*1, 56C 

BOX (I) *NB0X (I) »  <(G»ATl**2+GPAT2**2)*2.8E-2>/(GPAU»XlVGPATZ*X2) 
135  CONTINUE 

CALL  ERROR (POS, BOX, PH ],VNULL,NA) 

AN6LE*PHI 

CALL  ESTAN(ANGLE,PI ) 

ALPHA (2)* ANGLE 

EST*ANGLE 

OEG-EST 

CALL  ANOEG (OEG, PI ) 

OUHH(l)«OEG 

PRINT  *,  *EST  *  ',0EG 

00  1*0  X«l,*9 
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CALL  XPAT  (NR, ANGLE,  THETA, GP AT, G PAT 4 , GPAT8, GPAT  1, !PAT  2 
1  ,R,X1,X2,  DSEED,EST,  Vt*JLL»PI,NC»STF,  A) 

00  138  N-l, 568 

B0X<N>»NB0X(N>M(GPATl*»2PGPAT2**2>**.*E-2)/<CPAri*<l»GPAr2»X2) 
t>8  CONTINUE 

CALL  ERROR! POS, BOX, PHI ,VNULL«NA) 

8N6LE>PHI 

CALL  EST*N(AN6LE,PI> 

ALPMA!X«»>  ANGLE 
MTA-arrt*ALPHi(i*2> 

((f*KUTl«ESr 

ITMST 

CAU  *M(G(OE!,PI> 

IMM|*!>*OEt 

'wtUIUMt, ANEAN,8MEAN,VARi,VAR2,SNC,DUHM,NB> 

l»WIN 

HMi  1  >| I  ifMf 

••AM  ••  •••««*■• -2  /*AR2-»2 

•  •  I  ' 

Ul  COW'  to, 

RttUPw 

8*0 

IMKfTIW  *•*.*•  «  ».  T««TA,GPAT,  GRATA,  GPATB,G>»Tl,  SPAT? 

1  *  —  f  ,Ff  •  .PQS,M0X,NA> 

C  TM»  ••»•»*  rwt  PHASE  ERROR 

REAL  CTAttt  •,*•*«■  -*»•>«•> 

ANGLC-B. • 

CALL  ESTAM(ftMfc,.l.»l 

CALL  SOPAT!  anil l.T«r*  »  ..  . •  •  » .  -.-AT*,  JPNTl ,GPAT 2) 

SUM1-GPAT1 
SUH2-GPAT2 
00  220  I--260, 2 *1 
ANGLE-FLOAT (II  '110.  • 

P0S(I*28l)*I 

CALL  ESTANCANGLE, PI  I 

CALL  S0PAT!ANGLE,TNETA,GPAI»*PAT».  .— •  *  -  ff» 

EST1-GPAT1 

EST2-GPAT2 

ETA(I*2fll)« (SUN2*E5Tl”SUNl*EST2>  *2 «*4  - 1 
220  CONTINUE 

00  230  1*1, SGI 

MBOX<II-<ETA(I*l)-ETA(I>>/2.METAII) 

230  CONTINUE 

00  2*0  I- 1, 5E1 
POS(X)«POS(Z)'lOQ.t 
240  CONTINUE 
RETURN 
ENO 

SUBROUTINE  ESTANCANGLE, PI) 

C  THIS  PROGRAM  CONVERTS  DEG  TO  RAO 

ANGLE* AN! LE*PIF180 
RETURN 
ENO 

SUBROUTINE  SOP  A  T ( AN  GL  E ,T  NET A , GP  AT , GPAT  A,GPATB»GPATl,GPAT2l 
C  TNIS  PROGRAM  CALCULATES  THE  SUN  ANO  DIFFERENCE  AE  PATTERNS 
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THETA«»NSLEa.I  21598AA95 
CALL  PATT (THETA,  GPAT) 

GPATA«CPIT 

rHETI«IMCU>.  •  1 15  98  449B 

CALL  PATf (THETI»GPAT> 

6PATB*GPAT 

SPAT1»SPATA*6PATB 

GPAT2«CPATA-GPATB 

RETURN 

CNO 

SUBROUTINE  PAf T(THETA ,GPAT) 

IP  (THETt  • EO.  6)  THEN 
GPAT-1.0 
.  ,  ELSE 

. ,  ...  call  szh:<theta,gpat> 

. ENOiP 
'  RETURN 
eNO 

SUBROUTINE  SINC (THETA »GP AT) 

C  TNIS  PROGRAM  PIKVS  GPAT 

GPAT*SINU3.«  8  J  «97»  SIN  (THETA  > )  /  (A3. 1 8389T *SIN  (THETA) ) 

RETURN 

ENO 

SUBROUTINE  XPAT(NR, ANGLE, THETA, GPAT, GPATA,6PATB,;pIU, GPAT2 
1  *R,Xi,X2, OSEED,EST, VNULL, PI ,NC,STF, A) 

C  THIS  PROGRAM  CALCULATES  THE  XI  ,  X2  SIGNALS 

REAL  R(NC)»STF 
A*i.O 
ANGLE-1. s 

CALL  ESTAN (ANGLE, PI) 

CALL  SOPAT(ANGlE,THETA,GPAT,GPATA,GPATB,GPA(TL,GPAT2> 

CALL  NOISE(NR,t,OS£ED) 

XI-fl.0 

00  151  1*1, 1BC 

Xl*((R(I)»STF*GPATl)»A*2.8E-<t)*Xi 
151  CONTINUE 

CALL  NOISE ( NR, R,DSEED) 

*2*0.0 

00  1««  1*1, IBS  • 

X2«((R(I|*STF*GPAT2)*A»2.»E-4)*X2  . 

168  CONTINUE 
ANGLE*EST 

CALL  SOPAT (ANGLE, THETA, GPAT, GPATA,GPATB,GP ATI, GPAT 21 

WNULL*X1*  6PAT2-*2*G  PAT 1 

RETURN 

ENO 

SUBROUTINE  NOISC(NR,R, OSEEO) 

C  THIS  PNOGRAN  GENERATES  A  NOISE  VECTOR 

REAL  R (HO) 

OOUBLE  PRECISION  OSEEO 
NR* IN* 

CALL  GGNNL  (D$EEO,<tR,R> 

RETURN 

ENO 

SUBROUTINE  ERR  JR ( PO  S, BOX , PHI , VNULL , NA ) 

REAL  POS(NA),SOX(NA) 
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IF  (VNULL  .6E.  B0X(96«)>  THEN 
PHI-MS  <161  > 

goto  see 
ENOIF 

DO  250  1*1.960 

IF  (VNULL  «LE.  80X( X>  t  THEN 

PHI»MS(I> 

GOTO  2G0 
ENDIF 

250  CONTINUE 
260  CONTINUE 
RETURN 
END 

SUBROUTINE  ANOEGtDEG.PI) 

C  THIS  PROGRAM  CONVERTS  RADIANS  TO  OEGRBES 

OEG-DEGUSB/PX 
RETURN 
END 

SUBROUTINE  HEVIR< AL  PNA  .ANEAN.BNEAN.  VARi,  V  AR2,  SNC*  D'JHN.NB) 
REAL  ALPHA <NB)*OUNN<N8> 

ANEAN-O-I 
BMEAN-0.0 
VAR1-0.0 
VAR2-I.0 
00  135  I-lf SNC 
AMEAN* ALPHA (I+EI+AMEAN 
VARl*  ALPHA (I*2)**2*VAR1 
BHEAH-OUNMC I) ♦ SHE AN 
VAR2-0UMN (I l**2»VAR2 
13S  CONTINUE 

AHEAN-ANEAN/SNS 

VAR1-VAR1/SNOAHEAN** 2 

BMEAN-BNEAN/SNC 

VAR2*VAR2/SNC-BHEAN**2 

RETURN 

ENO 

SUBROUTINE  LOO ( Y, OH EA N,NB> 

REAL  V(NB).OMEAN(NB? 

00  270  1*1.50 

IF  (DNEAN(I)  .LE.  -3.S)  THEN 

OMEAN(X)—  3.0 

ELSE 

OMEAMdl-OHEANIX) 

ENOIF 

IF  (DHEAN(I)  .GE.  3.0)  THEN 
0MEAN(X)-3.« 

ELSE 

OMEAN(X)-OMEAN(I) 

ENOIF 

V  ( D-OMEAN  ( I) 

270  CONTINUE 
RETURN 
END 

SUBROUTINE  LOP< Y.BVAR  *NB) 

REAL  Y(NB)»8V1R(NB) 

DO  200  1*1.50 
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XF  (BVARC I)  .CE.  3.l>  THEN 

BVAR(X>>3.1 

ELSE 

BVAR(X)bBVAR(XI 

ENOXF 

r(2)*BVAt(I) 

260  CONTINUE 
RETURN 
END 

SUBROUTINE  LOOT r*MHU*NB> 
REAL  T(NB),MHU(NB) 

DO  All  1*1* 59 

XF  (MHU(I)  .GE.  3.1)  THEN 

MHU(X) *3.1 

ELSE 

NHU(X)bMHUTX) 

ENOXF 

V( X) *HHU( I) 

Alt  CONTINUE 
RETURN 
END 
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•  * 

•  PROGRAM  VELOEN  * 

•  » 

•  THIS  PROGRAM  PLOTS  THE  PROBALXTY  » 

•  DEMSZTr  FUNCTION  OF  THE  ANGLE  OFF  • 

•  BORESXGHT  FOR  A  MOVING  TARGET  » 

•  • 

•»**•»**»**»#**» •*•*••»•*••***•***••••••••••* 

PROGRAM  MAIN 

C  THIS  PROGRAM  SIMaATES  THE  SUFFICIENT  STASTI0 

INTEGER  NR,NA,MB,NC 

REAL  ANGLE, THETA, GPAT,GPATA,GPATB,GPAT1,GPAT2,PI,PHI 

1  , DEG, BE T A, AH E AN, B ME  IN, VARt , VAR2,SNC ,STF , OCC (SO, 311 , SORT (33) 

2  , EST,VNULL,SOHl ,SUH !,EST1, EST2,A ,HBOX (S61) ,MHU( 55)  , VEL 

3  , ETA(561> , POS <561) , BOX (561) , DIMM (5 3) ,ALPHAC5S» 

4  ,  R(100)  ,  XL,  X  2,  ONE  AN  (53)  ,  B  VAR  (53)  ,  XI 33) ,  V  (33)  ,ORS 

5  , EHEAN,CVAR 

6  ,  SNR 

DOUBLEPRECISION  OSEEO 

DATA  ANGLE,THET A, GP AT,GPAT A,GPATB»GPAT1,GPAT2,EST,PI 

1  , PHI, VNULL,SIIM1,SUMI,EST1,EST2, DEG, BETA, RhEAN.BHEAN, VAR1, VAR2 

2  , SNC,STF,X1, X  2, A, ORG.VEL 

3  /R,0,0,t,8,D,0,0«0«8,l,0,t,0,C,0,0,0,l,8,l,G,t,0,l,0,0,)f 
CALLPLOTS(0«,F.,9> 

CALL  PLOT<0.,-il.,-3> 

CALL  PLOT  <2«,S,5,-3) 

NA*561 

NB-53 

N5»100 

00  100  1*1, 1BG 
R(X)*0.0 
100  CONTINUE 

00  108  1*1,50 
DO  105  <1*1, 30 
OCCCI,J»*0.0 
105  CONTINUE 
108  CONTINUE 

00  110  1*1,53 
ALPHA (X) *0«  0 
OUMMm*I.O 
DHEAN(I>*0*0 
BVAR(l)*I.O 
MHU(I) *0,0 
110  CONTINUE 

00  115  1*1,33 

SORT(I)*0,0 

vm*o.o 

V(I)*0*0 
115  CONTINUE 

00  120  1*1,561 
ETA (I)*0,0 
POS(I) *0,0 
BOXIIMOiO 
MBOX(I) *0*0 
120  CONTINUE 
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DSEED»im57.D» 

CALL  FNECPI*ANCLE,THE1A,GPAT*GPATA,GPATB,  GPAT1,G*AT2 
1  *SUMlf SUM2*EST1,E  ST2*ETA  *POS*  MBOX ,NA) 
oo  ist  K«i,n 
ESTbO.B 
BETA-G.I 
STFaf, 

sNR>(LOGLi(iei.o/srF*»2.i))*io.e 
OSEED«DSEED*FL  OATCX) 

DO  125  I»3t|3liZ 

SORT ( CI+323/2) *FLOA  TCI)/ 10,0 

CONTINUE 

00  128  1*1* 30 

xmasoRTm 

SORT (I) a( SORT (I ♦1I~S0  *T( II ) 72* •♦SORT ( I ) 

CONTINUE 

CALL  ERIQCNR, ANGLE* THETA, CPAT*6PATA*GPATB*GPAT1, SPITZ* VEL 

1  ,R,X1*K2, DSEEO.EST, \NUU,PI, POS , BOX ,PHX ,  ALPHA* ARE AN* BHE AN 

2  ,VAR1,FAR2,DUNN,SNC  ,STF, NA ,NB ,NC , BETA, ONEAN ,BWAR, 980 X, MHO) 

00  LAO  1*1*50 

IF  CDUMMCI)  ,GE.  SORT (30 ))  THEN 
OCCCI,3B)«DCCCI,30) *1.8 
GOTO  A30 
ENOIF 

DO  A20  J*l, 30 

IF  CDUNHCI)  .LE.  SORT(J))  THEN 

OCC(I,JI*OCC(I*J) +1.8 

GOTO  A 38 

ENOIF 

CONTINUE 

CONTINUE 

CONTINUE 

CONTINUE 

DO  A8Q  L*l* 50 

ENEANal.l 

CVAR*0«C 

00  ATO  N«l* 30 

EMEAN*0C2 (L ,H) F50,3 ’SORT CH) +EMEAN 

CVAR*OCC ( L * N)  FIO.O*  CSORT CM) **2. I)  *CVAR  * 

CONTINUE  . 

CVAR«CVAR-ENEAN**2.0 

PRINT*  C/FT5 *A* T  28,F5, 2*T33*A//f T13, A*T29* I 2// 

1  T15,A,F10.2/FT11, A,F1C.2///T2,A/ 

2  T9,A/T8,A,T38,F5. 2,TA1 ,A/T 9, A ,T21 ,1 2///) *, 

3  'SIGNAL  TO  NOISE  RATIO  • », SNR, »0B », *AT  PULSE*, 

A  CL) , *HE AN  **,EHEAN, 'VARIANCE  «*,CVAR, 

5  *FIG  PLOT  OF  THE  PROBALITV  DENSITY  FUNCTION*, 

6  *0F  THETA  ESTIMATE  VERS  THETA*, 

T  'FOR  SIGNAL  TO  NOISE  RATIO  OF *,SNR, *08* , 

6  * AND  AT  PULSE** CL) 

X(31)«-3, 

YC31) «B  » 

XC32)*£.F6. 

VC32)*1.FA. 

00  A55  J*l,  30 
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Y( J)*0CC<1, JJ/IC.O 
W5  CONTINUE 

CALL  MIS  <«.,!. ,’THETA*,>8,6.,B.,X(3l)  ,X(32>> 

CALL  AXIS(3.,I.,*P(THETA) •,«8,A.,B8*,Y(3U  ,Y(32»I 
CALL  PLOT(8.,0.,3> 

CALL  LINE(X(V,38, 1,8,8) 

CALL  PLOT (8,, 8, ,3) 

CALL  PLOT(8,,-ll,,-3) 

CALL  PLOT(8.8«I«5,<*1) 

00  *80  1*18,88,18 
00  4$a  J-l, 38 

Y(J>*OCC(I,.»/fB.O 
A58  CONTINUE 

CALL  AXIS(0*,t •  »*THETA*,-5,6«»8»,X(31J*X(32)) 

CALL  AXIS  (3 *,8  »,*P( THETA) '  ,*6  ,*  •  ,90«,  T  (31) ,T(32>> 

CALL  PLOT (8  •,!• ,3) 

CALL  L1NE(X,Y, 38, 1,0,0) 

CALL  PLOT(8*,8.,3> 

CALL  PLOT  (0., -11., -3) 

CALL  PLOT  <8.0*1  .8,-3) 

*69  CONTINUE 
*80  CONTINUE 

CALL  PLOTE(N) 

ENO 

SUBROUTINE  ERI 0 (NR, ANGLE, THETA, GPAT ,GPATA , GPATB, SPAT  1 ,  GPAT2 ,V EL 

1  ,R«Xl,K2, DSEEO, EST, VNULL,PZ,POS, BOX , PHI, ALPHA, A MEAN, BHEIN 

2  ,VAR1,VARZ,DUMM,SNC,STF, NA,NB »NC, BET  A,DMEAN  *BVAR,HB0X ,MHU) 
REAL  POS(NA),80X(NA), CUNN(NB), ALPHA (NB),HHU(NB> 

1  , R (NCI , OHEAN! NB) , BVARI NB) , NBOX (NA  > 

VEL*1.8 

CALL  XPAT (NR, ANGLE, THETA ,6PAT »6PAT  A,6PATB» GPAT 1, SPAT  2 
L  ,R,X1,K2,0S£ED,EST, VNULL»PI»NC,STR» A»VEL> 

00  135  I«1,56C 

B0X(I>*HB0X(I>»((GPATl*»2*GPATZ**Z>*2.9E-2>/(GPAU»Xl»GPAT2*X2> 
135  CONTINUE 

CALL  ERROR! POS , BOX, PHI , VNULL ,NA ) 

ANGLE*PHI 

CALL  ESTA N (ANGLE, PI > 

ALPHA !2»>ANGLE 

EST*ANGLE 

0E6*EST 

CALL  ANDEG(OEG,PI) 

OUHNI1JOEG 
DO  1*9  1*1, 49 
VEL*FLOAT(I)U.C 

CALL  XPAT (NR, ANGLE, TMETA, GPAT,6PATA,GPATB,GPAT1,5PAT2 
1  ,R,X1,X2, DSEEOyEST, VNULL ,PI,NC,STF, A ,VEL ) 

DO  138  N*l,560 

BOX (N) *N80X (N) * ( (GP AT 1*» 2*GPAT2**2>*2.8E-2 )/ ( GPATl* X1*GPAT2»X2> 
138  CONTINUE 

CALL  ERRDR( POS, BOX, PHI, VNULL ,NA) 

ANGLE-PHI 

CALL  ESTAN(ANGLE,PII 
ALPHA  «♦« -ANGLE 
EST*EST ♦ANGLE 
OEG-EST 


CALL  ANDES (DEG, PI) 

0UMH(I41)*0EG 

SNC«Z4& 

CALL  HEVAR(ALPNA , ANEAN,BMEAN, VARi, VAR 2, SMC ,0UMM,V9) 
0HEAN(I*O«BMEAN 
B VAR (141) -VAR2 
HHU  C  X  41 )  «  BNEAN  »  *  2 /V  AR  2**  2 
14#  CONTINUE 
RETURN 
ENO 

SUBROUTINE  FNE  ( PI ,A NGLE,THETA,GPAT, GPATA,GPATB,G*AT1,CPAT2 
1  ,SUN1,SUH2,EST1 ,EST J »ETA ,POS,HBOX ,NA) 

C  THIS  PROGRAM  ESTIMATES  THE  PHASE  ERROR 

REAL  £TA(NA),P9S(NA ),HBOX(NA) 

ANGLE*!. I 

CALL  ESTAN(ANGLE,PX ) 

CALL  SOPAT( AN«LE,THETA,SPAT,GPATA,GPATB,SPffTl,GPlT2) 

SUH1-6PAT1 

SUN2-GPAT2 

00  221  I— 281,201 

ANGLE-FLOAT (X) /ISO. 0 

POS (14281) -I 

CALL  ESTANIANGLE, PX ) 

CALL  SOPAT (ANGLE, THETA, SPAT, GPATA,GPATB,GPATi,GPlT2l 

EST1-GPAT1 

EST2-GPAT2 

ET A (I42S1 1  *(SU*2*E$Tt •SOMl*E$T2>*2*OE-2 
220  CONTINUE 

DO  230  1-1,56# 

NBOX  (I)-(ETA(I4U-ETA(t))/2.P+ETA(I) 

230  CONTINUE 

00  240  1-1,561 
POS(I)-POS(I)/lOO.O 
240  CONTINUE 
RETURN 
ENO 

SUBROUTINE  ESTAN( ANGLE, PI) 

C  THIS  PROGRAM  CONVERTS  DEG  TO  RAD 

angle«ansle*pi/ibo 

RETURN 

ENO 

SUBROUTINE  SOPATIANGL E, THETA, GPAT,GPATA,GPATB,GPATt,GPAT2) 

C  THIS  PROGRAM  CALCULATES  THE  SUM  AND  DIFFERENCE  AE  PATTERNS 

THETA-ANSL £4,1215904495 
CALL  PATT (THET  A , GP  AT ) 

GPATA-GPAT 

THETA- ANGLE'* 1 2 15 98 4495 
CALL  PATT (T)CTA,GPA  T) 

GPATB-CPAT 

SPAT 1-SPAT A  4CPATB 

6PAT2-GPATA-6PATB 

RETURN 

ENO 

SUBROUTINE  PATT (THETA ,GPAT) 

IF  (THETA  ,E0.  B)  THEN 
CPAT-1,0 
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CALL  S1N3 (THETA, GPAT) 

ENOIF 

RETURN 

END 

SUBROUTINE  SINS (THETA, GPAT) 

THIS  PROGRAM  FINDS  GPAT 

GPAT>SIM(R3«0P  3097*  SIN  (THETA) )  /  (N3.Q03O9E  *  SIN  (THET  A)  > 

RETURN 

ENO 

SUBROUTINE  XPAT (NR, ANCLEfTHETA, GPAT , GRATA, GPATB, GRATA, GPAT 2 
1  ,R,K1,K2,OSEEO,EST, \NULL ,PI,NC,STF, A, VEL) 

THIS  PROGRAM  CALCULATES  THE  XI  ,  XZ  SIGNALS 
REAL  R(N3) ,STF 
A«1.8 

ANGLE* VEL /1 0.1 -Z.  0 
CALL  ESTAN(ANGLE,PI) 

CALL  S0P1T(ANGLE,THETA,GPAT,GPATA,GPATB,GPAT1,GP1T2> 
call  NOISE(NR,R,DSEEO> 

Xl«C.O 

oo  tsa  i«i,iog 

Xt»( (R(I)  *STF*GPAT1 )*A*2»0E“4) *X1 
CONTINUE 

CALL  NOI$E(NR,R,DSEED) 

X2«C.O 

00  160  1-1,100 

X2«((R<II*STF*GPAT2>*A*2.0t-A)*tt 

CONTINUE 

ANGLE*EST 

CALL  SOPAT (ANGLE, THETA, GPAT, GPATA,GPATB,GPATl,GPAt2t 

VNULL*X1*GPATZ-XZ*GPAT1 

RETURN 

ENO 

SUBROUTINE  NOISE (NR ,R  ,DSEEO) 

THIS  PROGRAM  GENERATES  A  NOISE  VECTOR 
REAL  R (100) 

00U6LE  PRECISION  OSEED  ' 

NR-189 

CALL  G6NML  (OSEEO,NR,R) 

RETURN 

END 

SUBROUTINE  ERROR ( POS, EOX ,PHI ,VNULL ,NA) 

REAL  POS(NA) ,BOX(NA ) 

IF  ( VNULL  ,GE.  BOX(  96  )))  THEN 
PHI-POS(f 51) 

GOTO  260 
ENOIF 

00  290  I*l«96l 

IF  (VNULL  «LEi  BOX ( I) )  THEN 

PHI-POS(I) 

GOTO  269 
ENOIF 
CONTINUE 
CONTINUE 


RETURN 

END 

SUBROUTINE  ANDEGIDEG,PI> 

THIS  PRO* RAN  CONVERTS  RADIANS  TO  DEGREES 
DES*DE6*t6B/PI 
RETURN 
ENO 

SUBROUTINE  HEVAR(ALPHA,ANEAN»BMEANf VAR1, VARE*  SNC» OUNH(N6) 
REAL  ALPRA(NB) »OUHH INB) 

ANEANM.I 

BNEANM.B 

VARL«a.B 

VAR2-0.0 

00  L3S  Z*1,SNC 

AHEAN*ALPHA II*E IPAHEAN 

VAR1*  ALPHA  II ♦VAR1 

BNEANaOUNN(Z)»BNEAN 

VAR2>0UM*im»»E*VAR2 

CONTINUE 

ANEAN*ANEAN/SNC 

VKR1«VARL/SNC*INEAN*>  2 

BNEAN*BNEAN/SNC 

VAR2»VAR2/SNC-§MEAN**2 

RETURN... 

ENO 
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